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Abstract: Within the Quantum Field Theory context the idea of a "cosmological con- 
stant" (CC) evolving with time looks quite natural as it just reflects the change of the 
vacuum energy with the typical energy of the universe. In the particular frame of Ref. 
1 31], a "running CC" at low energies may arise from generic quantum effects near the 
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Planck scale, Mp, provided there is a smooth decoupling of all massive particles below 
Mp. In this work we further develop the cosmological consequences of a "running CC" by 
addressing the accelerated evolution of the universe within that model. The rate of change 
of the CC stays slow, without fine-tuning, and is comparable to H 2 Mp. It can be described 
by a single parameter, u, that can be determined from already planned experiments using 
SNe la at high z. The range of allowed values for v follows mainly from nucleosynthesis 
restrictions. Present samples of SNe la can not yet distinguish between a "constant" CC 
or a "running" one. The numerical simulations presented in this work show that SNAP 
can probe the predicted variation of the CC either ruling out this idea or confirming the 
evolution hereafter expected. 
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1. Introduction 



The Standard Cosmological Model fits our universe, in the large, into an homogeneous 
and isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmological type [jl]]. Its 
4-curvature is determined from the various contributions to its total energy density, namely 
in the form of matter, radiation and cosmological constant. 

Evidence for a dominant content of energy in the form of cosmological constant was 
found by tracing the rate of expansion of the universe along z with high-z Type la su- 
pernovae [§, |^|. This measurement combined with the measurements of the total energy 
density Qj, from the CMB anisotropies [§], H, indicates that ~ 70% of the critical energy 
density of the universe is cosmological constant (CC) or a dark energy candidate with a 
similar dynamical impact on the evolution of the expansion of the universe. The matter 
content, on the other hand is dominated by the dark matter, whose existence is detected 
by dynamical means |T|, and amounts to Q° M ~ 30% of the critical density. 

The CC value found from Type la supernovae at high z & || is: 

A = 0^p°~ 6 hi xHT 47 GeV 4 . (1.1) 

Here p° c ~ (3.0 V^-o x 10~ 12 GeV) 4 is the present value of the critical density, and ho ~ 0.7± 
0.1 sets the typical range for today's value of Hubble's constant Hq = 100 (Km/ sec Mpc) ho. 
In the context of the Standard Model (SM) of electroweak interactions, this measured CC 
should be the sum of the original vacuum CC in Einstein's equations, A mc , and the in- 
duced contribution from the vacuum expectation value of the Higgs effective potential, 
A«nd = (V e g): 

A = A vac + A ind . (1.2) 

It is only this combined parameter that makes physical sense, whereas both A vac and Ai n< i 
remain individually unobservable 1 . From the current LEP 200 numerical bound on the 
Higgs boson mass, M H > 114.1 GeV §, one finds \A ind \ = (1/8)M^ v 2 > 1.0 x 10 8 Gel/ 4 , 
where v ~ 246 GeV is the vacuum expectation value of the Higgs field. Clearly, \A inc i\ is 
55 orders of magnitude larger than the observed CC value ( |1.1| ). Such discrepancy, the 
so-called "old" cosmological constant problem |?| ||, manifests itself in the necessity of 
enforcing an unnaturally exact fine tuning of the original cosmological term A vac in the 
vacuum action that has to cancel the induced counterpart Aj„^ within a precision (in the 
SM) of one part in 10 55 . 

The measured CC remains very small as compared to the huge CC value predicted 
in the SM of Particle Physics. Actually, if the physical value of the CC would conform 
with that one predicted in the SM, the curvature of our universe would be so high that the 
Special Theory of Relativity could not be a solution to Einstein equations to any reasonable 
degree of approximation. Therefore, the SM prediction of the CC violently contradicts our 
experience, whereas the small measured value Ql.l|) is perfectly compatible with it. 



1 In general the induced term may also get contributions from strong interactions, the so-called quark 



and gluon vacuum condensates. These are also huge as compared to (1.1), but are much smaller than the 
electroweak contribution V e a. 
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The Cosmological Constant Problem (CCP) is a fundamental problem. It is most 
likely related to the delicate interplay between Gravity and Particle Physics, and it has 
become one of the main poles of attention ||. All attempts to deduce the small value 
of the cosmological constant from a sound theoretical idea ended up with the necessity 
of introducing severe fine-tuning. This concerns also, unfortunately, the use of supersym- 
metry and string theory (see e.g. H In this respect we recall that, for a realistic 

implementation of the existing versions of M-Theory, one would like to have a negative 
(or at least vanishing) cosmological constant in the remote future, such that it does not 
prevent the construction of the asymptotic S-matrix states in accelerated universes [11]. 
Since the presently observed value of the CC is positive, there is the hope that a variable 
cosmological term may solve this problem. 

There is a permanently growing flux of proposals concerning the CCP. On the first place 
there is the longstanding idea of identifying the dark energy component with a dynamical 
scalar field [12, 13]. More recently this approach took the popular form of a "quintessence" 
field slow-rolling down its potential [14] . This proposal has, on its own, given rise to 
a wide variety of models JO], |lq] . Extended models of this kind ( "k-essence" ) are also 
based on scalar fields but with a non-canonical kinetic energy [17|. The main advantage 
of the quintessence models is that they could explain the possibility of a variable vacuum 
energy. This may become important in case such variation will be someday detected in 
the observations. Recently other approaches have appeared in which the dark energy is 
mimicked by new gravitational physics |l8|]. From the point of view of the CCP, all these 
approaches lead to the introduction of either a very small parameter or a very high degree 
of fine-tuning. In another, very different, vein the possibility to accept the observed value 
of the CC within the context of a many world pool is offered by the anthropic proposal |l9|| . 
Let us finally mention the intriguing proposal of non-point-like gravitons at sub-millimeter 
distances suggested in p0| , or the possibility of having multiply degenerate vacua |2lf| . 

When assessing the possibility to have variable dark energy, other no less respectable 
possibilities should be taken into account. In a series of recent papers ]22|, [23]] , the idea 
has been put forward that already in standard Quantum Field Theory (QFT) one would 
not expect the CC to be constant, because the Renormalization Group (RG) effects may 
shift away the prescribed value, in particular if the latter is assumed to be zero. Thus, 
in the RG approach one takes a point of view very different from e.g. the quintessence 
proposal, as we deal all the time with a "true" cosmological term. It is however a variable 
one, and therefore a time-evolving, or redshift dependent: A = A(z). Although we do not 
have a QFT of gravity where the running of the gravitational and cosmological constants 
could ultimately be substantiated, a semiclassical description within the well established 
formalism of QFT in curved space-time (see e.g. [24, ^]) should be a good starting point. 
From the RG point of view, the CC becomes a scaling parameter whose value should be 
sensitive to the entire energy history of the universe - in a manner not essentially different 
to, say, the electromagnetic coupling constant. One of the main distinctions between our 
approach and all kinds of quintessence models is that these models imply the introduction 
of a phenomenological equation of state p x = w p x for the scalar field x mimicking the CC, 
where w is a negative index (smaller than —1/3). Whether constant or variable, a "true" 
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cosmological parameter has, instead, no other equation of state associated to it apart from 
the exact w = — 1 one. 

Attempts to apply the RG for solving the CC problem have been made earlier [26, 
27]. The canonical form of renormalization group equation (RGE) for the A term at high 



energy is well known - see e.g. [25, |28|. However, at low energy decoupling effects of 
the massive particles may change significantly the structure of this RGE, with important 
phenomenological consequences. This idea has been retaken recently by several authors 



from various interesting points of view [22, ^j, |2^, |3(J. However, it is not easy to achieve 



a RG model where the CC runs smoothly without fine tuning at the present epoch. In 



Ref. [31, |32[ a successful attempt in this direction has been made, which is based on possible 
quantum effects near the Planck scale. At the same time, the approximate coincidence of 
the observed A and the matter density, ~ ^j}/ ' 1 - e - the " new " CC problem, or "time 
coincidence problem" || can be alleviated in this framework if we assume the standard 
(viz. Appelquist-Carazzone [^]) form of the low-energy decoupling for the massive fields. 

In the present paper we elaborate on this idea further. We develop a semiclassical 
FLRW model whose running CC is driven smoothly, without fine tuning, due to generic 
quantum effects near the Planck scale. We show that, due to the decoupling phenomenon, 
the low-energy effects (in particular the physics from the SM scale) are irrelevant for the CC 
running, and so the approximate coincidence between £l® M and is not tied to any par- 
ticular epoch in the history of the universe. Furthermore, the new effects imply deviations 
from the standard cosmological equations due to quantum effects. Our "renormalized" 
FLRW model provides a testable framework that can be thoroughly checked from SNAP 



data on Type la supernovae at very high z - see [34, 35] and references in 36]. If these 
experiments detect a z-dependence of the CC similar to that predicted in our work, we may 
suspect that some relevant physics is going on just below the Planck scale. Alternatively, 
if they find a static CC, this might imply the existence of a desert in the particle spectrum 
near the Planck scale. 

The structure of the paper is as follows. In the next section we compare constant versus 
variable CC models. In Section |3] we present our variable CC model based on the Renorma- 
lization Group. In Section || we solve the FLRW cosmologies with running CC. In Section 
|5| we study the numerical behaviour of these cosmologies, and the predicted deviations 
from the standard FLRW expectations. In Section ^, we introduce the magnitude-redshift 
relation for the analysis of the SNe la. In Section ^ we perform the simulations on the 
SNAP data in order to test the sensitivity with which the features of the new model can 
be determined. In the last section we draw our conclusions. Two appendices are included 
at the end: one to discuss some technical issues inherent to our QFT framework, the other 
providing some background related to the statistical analysis. 

2. Constant versus variable cosmological term 

The cosmological constant enters the Hilbert-Einstein (HE) action as follows 

She = ~ I d A x^r~g (j^ R + A vac) ■ (2.1) 
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It is well-known that renormalizability requires that this effective action should be extended 
with a number of higher derivative terms [24, 25 1: 

S vac = J d A Xsf^) { aiRl vap + a 2 Rl u + a 3 R 2 + a A BR - ^ Q R - A mc } . (2.2) 

The phenomenological impact of the higher derivative terms in this action is negligible at 
present, and therefore it suffices to use the low-energy action ( |2.1| ). However, the presence 
of the parameter A vac is as necessary as any one of these higher derivative terms to achieve 
a renormalizable QFT in curved space-time 2 . 

The vacuum CC itself, A vac , is not the physical (observable) value of the cosmologi- 
cal constant. By definition the physical CC is the parameter A entering the Friedmann- 
Lemaitre equation: 



d\ 8-7T Gn k 



where H is the expansion parameter (Hubble's "constant" ) . This A should be the sum ( |l,2j ) . 
While the homogeneous and isotropic FLRW cosmologies do not allow spatial gradients of 
A, they do not forbid the possibility that A may be a function of the cosmological time: 
A = A(i). In this case the Einstein field equations associated to the action ( |2.1| ) read 



R»v - \g^R = -8irG N , (2.4) 

where T^ u is given by T^ v = T^ u + g^ u k{t), T^ u being the ordinary energy-momentum 
tensor associated to matter and radiation. By the Bianchi identities, it follows that A is 
a constant if and only if the ordinary energy-momentum tensor is individually conserved 
(S/^T^y = 0). In particular, A must be a constant if T^ u is zero (e.g. during inflation). 

Modeling the expanding universe as a perfect fluid with velocity 4- vector field U^, we 
have 

T^ = - P g^ + (p + p)U»U\ (2.5) 

where p is the isotropic pressure and p is the proper energy density of matter. Clearly 
takes the same form as ( p.5| ) with p — > p = p + A , p ^> p = p — A. Using the FLRW 
metric 

ds 2 = dt 2 - a 2 (t) ( dr i = + r 2 d6 2 + r 2 sin 2 9 d<p 2 ) , (2.6) 
\ 1 — k r z J 

we can compute explicitly the local energy-conservation law Tn V = 0. The result is the 
old Bronstein's equation [37] allowing transfer of energy between ordinary matter and the 



dark energy associated to the A term: 

A + p + 3H(p + p) = 0. (2.7) 

We see that the most general local energy conservation law (or equation of continuity) 
involves both the time evolution of p and that of A. For a truly constant CC, then A = 0, 



2 It follows that quintessence models without a A term cannot be renormalizable theories in curved 
space-time. 
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and we recover of course the standard conservation law p+3H (p+p) = 0. Equations (|2.3| ) 
and ( |2.7D constitute two independent counterparts for constructing FLRW cosmologies with 
variable A. The dynamical equation for the scale factor is 

47T 47T 

5= — -G N (p + 3p)a = — -G N (p + 3p-2A)a, (2.8) 



but it is not independent from the previous two. In the matter era p = 0, and Eq. 
shows that unless |A| is much smaller than p, a positive A eventually implies accelerated 
expansion - as in fact seems to be the case for our universe |2], [|. 

It should be clear that our approach based on a variable CC departs from all kind of 
quintessence-like approaches, in which some slow-rolling scalar field \ substitutes for the 
CC. In these models, the dark energy is tied to the dynamics of x, whose phenomenological 
equation of state is defined by p x = w p x . The term —2 A on the r.h.s. of Eq. (|2.8j) must 
be replaced by p x + 3p x = (1 + 3u>)p x . In order to get accelerated expansion in an epoch 
characterized by p = and p —* in the future, we must require — u>_ < w < —1/3, where 
usually u>__ > — 1 in order to have a canonical kinetic term for \- However, one cannot 
completely exclude "phantom matter-energy" {w- < — 1) and generalizations thereof p8[ . 
Present data suggest the interval —1.38 < w < —0.82 at 95% C.L. [39]. Although p x and 



p x are related to the energy-momentum tensor of x, the dynamics of this field is unknown 
because the quintessence models do not have an explanation for the value of the CC. 
Therefore, the barotropic index w is not known from first principles. In particular, one 
cannot exclude it may have a redshift dependence, which can be parametrized in various 
ways as follows: 

E2L = w = Wo + WlZ + o(z 2 ) = w + w a (l - a) + 0((1 - a) 2 ) , (2.9) 
Px 

where z + 1 = 1/a. Finding a non- vanishing value of w\ (w a ) implies a redshift evolution 
of the equation of state for the x field || . The difficulties to measure w\ are well-known, 
see e.g. fitf . 

Quite in contrast to that scenario, since our variable CC is a "true" cosmological 
parameter, the only possible equation of state for the CC term is p = —A, whether it is a 
true constant or it is a parameter that evolves with the cosmological time. In our case the 
CC is indeed a variable one, and its variation is attributed to potential quantum effects 
linked to physics near the Planck scale, as will be explained in the next section. 

3. Renormalization group and cosmological constant 

The possibility of a cosmological model with a time-dependent A as presented in the previ- 



ous section is very generic. However, the two differential equations ( |2.3| ) and (2.7) cannot 
be solved unless a third equation involving A is called for. The third equation admits many 
formulations, even at the classical level 3 . However, a particularly interesting implemen- 
tation occurs when the time dependence has its prime origin in the quantum field theory 



3 See e.g. ^] and references therein. 
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notion of Renormalization Group running 23 j. It means that the Hilbert-Einstein 
action ( |2.1| ) is treated semiclassically and one introduces an equation for the running cos- 
mological constant. Although this can be done in several ways, a consistent formulation of 



the approach has been presented in [22] within the well established formalism of QFT in 



curved space-time (see e.g. [24, pq]). From simple dimensional analysis, and also from dy- 



namical features to be discussed below, the RGE for the physical CC may take in principle 



the generic form [22 



^ S>-f + /^>«? W£§ + ... 

\ * 3 3 3 3 

oo / \2n 

= EE«-^i ^ =/3 A (M,M/M). (3.1) 

where the sums are taken over all massive fields; A, B,C, D, ... are constant coefficients, 
and fi is the energy scale associated to the RG running. We assume that fi is of the order 
of some physical energy-momentum scale characteristic of the cosmological processes, and 
can be specified in different ways (see below). In our model we assume that [x is given by 
the typical energy-momentum of the cosmological gravitons, namely \i = H, which is of 
order R 1 / 2 . The r.h.s. of ( |3.1| ) defines the /^-function for A, which is a function of the 
masses and in general also of the ratios of the RG scale and the masses. 

In the equation above the masses of the various degrees of freedom (d.o.f.) are repre- 
sented by rrii and Mj. Here we distinguish between the active (or "light") d.o.f. at the scale 
H, namely those satisfying \x 3> m, and contributing to the /^A-function in the form ~ mf, 
from the "decoupled" (or "heavy") d.o.f. which satisfy fi <C Mj and yield the remaining 
terms in the series expansion of f3\ in powers of n/Mj <C 1. As can be seen, all the terms 
in Pa are of the form \j? n M.\~ 2n (n = 0, 1, 2, 3...) where Mi = mj for n = and Mi = Mi 
for n > 1. The n = terms correspond precisely to the active d.o.f. contributing the full 
fourth power of their masses. The coefficients Ai for these terms are known in the ultra- 
violet (UV) regime because they must coincide, in any mass-dependent renormalization 
framework, with their values in the Minimal Subtraction (MS) scheme |43| , fHf . In this 
regime the /^-function depends only on the masses of the active degrees of freedom, and 
not on the ratios ///Mj of the RG scale and the heavy masses. For particles i = 1,2,3, ... 
of masses mi and spins Ji one finds [22]: 



Ai = (-l) 2Ji (Ji + 1/2) nji N c , (3.2) 

with n{ ,i,i/2} = (1) 1) 2) and N c = 1, 3 for uncolored and colored particles respectively. The 
remaining terms in ( |3.l| ) "decouple" progressively faster as we move from n = 1,2,3, ... 

Notice that dimensional analysis is not enough to explain the most general structure 
of f3\. The fact that only even powers of /i are involved stems from the covariance of the 
effective action. Indeed, the odd-powers of \x cannot appear after integrating out the higher 
derivative terms, as they must appear bilinearly in the contractions with the metric tensor. 
In particular, covariance forbids the terms of first order in fj,. As a result the expansion must 
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start at the \i -order. On the other hand, the structure for the n > 1 terms associated to the 



coefficients B, C, ... in (3.1) is dictated by the the Appelquist-Carazzone (AC) decoupling 
theorem [33, |43|] . Thus, when applying the AC theorem in its very standard form to the 
computation of (3\, the decoupling does still introduce inverse power suppression by the 
heavy masses (those satisfying Mj 3> fj,), but since the /?A-function itself is proportional 
to the fourth power of these masses it eventually entails a decoupling law l/Mj" --4 , and 
so the n = 1 and n = 2 terms do not decouple in the ordinary sense whereas the n > 3 
terms do, i.e. only the latter start getting (increasingly higher) inverse power suppression 
by the heavy masses. The upshot is that, strictly speaking, the truly decoupling terms 
in f3\ (in the sense used when applying the AC theorem to the ordinary SM interactions) 
commence at n = 3 and above. In contrast, the n = 2 terms are constant (independent of 
the masses) and the n = 1 terms acquire the peculiar structure fj 2 M? hence displaying 
the unusual property that a /^-function may increase quadratically with the heavy masses 
("soft decoupling"). Remarkably enough, the CC is the only parameter in the effective 
action of vacuum that has the necessary dimension to possess this distinctive property 



22, pSfl , and the latter is certainly not shared by any other parameter in the SM of the 
strong and electroweak interactions. 

Despite that the explicit derivation of the decoupling for the CC is not possible at 
present (see the extended discussion of this issue in the Appendix 1), the assumed form of 



decoupling is highly plausible [22, |2J| within the general effective field theory approach ff5[ . 
Furthermore, the /^-functions for the remaining coefficients at (k = 1,2,3,4) of the vac- 
uum effective action (viz. those corresponding to the higher derivative terms in Eq. Q2.2j)) 



do exhibit exactly this kind of decoupling behaviour assumed for /?a [ 46 1 . In this situation 
it is quite reasonable to apply the phenomenological approach. Since there are no direct 
theoretical reasons to exclude the soft decoupling in the CC sector, we just admit that it 
really takes place and investigate the cosmological model which follows from this assump- 
tion. We will see indeed that the n = 1 structure (undoubtedly the most peculiar one 
of the /?A-function) can be experimentally probed in the next generation of high redshift 
cosmological measurements |34j]- Cf. Sections 

Following the phenomenological indications, it is very important that the structure of 
the /?A-function does not trigger a too fast running of A, which would be incompatible with 



the present observations 0, H]. From Eq. (3.1) it is clear that this feature will depend not 



only on the values of the masses of the various d.o.f. involved, but also on the characteristic 
energy scale \x used to track the RG running, which must be correctly identified. This is 
particularly evident from the quadratic structure of the n = 1 terms. In the following we 
mention a few different scenarios that have been contemplated in the literature: 

• In Ref. [^] it was assumed that only the lightest d.o.f. would contribute, equivalently 
Bj = Cj = Dj = ... = in Eq. ( |3.1| ). The only non- vanishing coefficients here are 
some of the Ai ^ 0, namely those associated to d.o.f. for which [i > rrii. Typically, 
this would be the case for the lightest neutrinos, whose mass can border the range 



m u ~ 10 eV [47], and therefore satisfy the curious numerical coincidence Aq ~ m 



4 



which motivated the RG approach of p3]. Moreover, in this paper the RG scale 
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was identified from the value of the fourth root of the critical density at a given 
cosmological time t: 

M~Pc /4 (*), (3-3) 

For the present universe, this scale is (p®) 1 ^ ~ W~ 3 eV, i.e. of the order of the 
lightest neutrino mass mentioned above. For the radiation era, p c ~ T 4 and so p in 
that epoch is essentially given by the temperature (p ~ T) within this Ansatz. 

In Ref. 1 26] the same RG scale ( |3.3| ) was adopted, but the important point was made, 
on the basis of effective field theory arguments, to the necessity of including the heavy 
mass terms Mj in Notwithstanding, when applying this framework to the SM of 
the strong and electroweak interactions, where the largest masses are of the order of 
a few hundred GeV, one is forced to tame the runaway evolution of A - triggered by 
the quadratic n = 1 terms on the r.h.s. of Eq. (|3.l|). In practice, it means that one 



has to enforce a fine tuning of their overall effect to zero [29| 

^ ±Jj (J, ivJ 

3 



J2B j ^Mf = (p°) 1/2 J2BjMj = Q, (3.4) 



otherwise one gets an extremely fast running of the CC which would be incompatible 
with the observations [Q, J3|. The authors of use this adjustment to hint at the 
value of the Higgs mass Mh, which (for particular values of the coefficients Bj in 
their given setting) is the only free mass parameter in the sum ( |3.4| ), that runs over 
all SM particles. The result that they obtain is reasonable (Mh ~ 550 GeV), but 
still too high as compared to the current bounds and expectations [^, ^8|. Moreover, 
the obtained value for Mjj is scheme-dependent. 

Eq. fl3,l| ) was proposed in Ref. [22] assuming that the RG scale is identified with the 
square root of the curvature scalar p ~ i? 1 / 2 , which in the FLRW cosmological context 
is equivalent to identify p with the expansion parameter (or "Hubble constant") at 
any given cosmological time: 

p ~ H(t) . (3.5) 

For the present universe, Hq ~ 10 -33 eV. This scale is much smaller than ( |3.3[ ), but 
from our point of view is the most natural one, as it is naturally linked with the scale 
of the cosmological gravitational quanta (gravitons) - used here in a generic sense 
referring to the presumed quanta of gravity as a field theory with a tensor potential, 
rather than to its relation with the gravitational waves. Scale (13.51) is also used 



successfully in other frameworks, e.g. in [49] to describe the decoupling of massive 
particles in anomaly-induced inflation. 

In Ref. Q] the RG scale p was identified with the inverse of the age of the universe 
at a given cosmological time, i.e. p ~ 1/t. This is essentially equivalent to the 
previous case, because H ~ 1/t in the FLRW cosmological setting. Nevertheless, 



the constitutive relations for the RG evolution in Ref. [30] are different from [22] 



and they are phrased in a non-perturbative quantum gravity framework based on the 



-8- 



(hypothetical) existence of an infrared (IR) fixed point. In our case, the RG approach 
aims at the simplest possible modification of the FLRW cosmology, namely the study 
of the CC evolution within perturbative QFT in a curved background. In contrast 
to [|30|] , we allow transfer of energy between matter /radiation and CC, but we do not 
consider any significant scaling evolution of Newton's constant. Indeed, in [22] it was 



shown that Gm does not undergo any appreciable running within our perturbative 
framework. 

• Finally, we consider the framework which we will elaborate in the rest of this paper. 
It is based on the identification (|3.5| ) and assumes that the heaviest possible masses 



entering Eq. ( |3.1[) lie near the Planck scale, Mp [31, 32]. This approach does not 



have any fine-tuning problem in the value of /3a, as we shall see. 



In the last framework the RGE that supplements fl2.3|) and (|2.7| ) is given by a particular 
form of Eq. ( |3.lD , namely 

= 7A2 E « +•••• = 7A2 E « r2m ? + ••• • ( 3 - 6 ) 

din /J, (47r) 2 (47r) 2 ^ 

where Mj is a collection of (superheavy) sub-Plan ckian-size masses just below the Planck 
scale, Mf < Mp. The remaining masses are the set of "low-energy" masses, mj, in the 
sense that mf <C Mf, and therefore do not contribute in any significant way to this RGE. 
Since fi = H is so small at present, there is not a single d.o.f. satisfying fj, > m, i.e. all 
coefficients A{ for the n = terms in Eq. ( |3,1| ) are zero. Then all the masses are supposed 
to "decouple" according to the soft terms 4 (~ H 2 Mf). 

Looking at the decoupling law, the n = 1 soft decoupling terms are always the leading 



ones as compared to all others on the r.h.s. of Eq. (3.1). Notice that we assume H < Mi and 
that the physics of Planckian or trans-Planckian energies is governed by some unspecified 
more fundamental framework (e.g. string/M theory). Therefore, we do not address here 
the issue of whether trans-Plankian physics may also be responsible for the CC or dark 
energy in general [ |50|| . In fact, we rather propose that the main contribution to the CC 
at present can be the sole result of quantum effects from the highest possible, but still 
sub-Planckian, energy scales. 
We also note that 

cHlMf/{A-K) 2 ~ c (1.5 10~ 42 Gey x 1.2 10 19 GeV) 2 / (47r) 2 ~ 10~ 47 GeV 4 ~ A (3.7) 

for some c = 0(1 — 10) and Mi ~ Mp. The previous result is very close to the observational 
data H ||]. This is highly remarkable, because two vastly different and (in principle) totally 
unrelated scales are involved to realize this "coincidence": Hq (the value of /i at present) 
and Mp, being separated by more than 60 orders of magnitude. In other words, the 
"coincidence" amounts to saying that the mass scale associated to the CC at present, 



4 Let us notice that the present-day Hubble parameter, Ho ~ 10 -33 eV, is 30 orders of magnitude smaller 
than the mass of the lightest neutrino, 41 orders of magnitude smaller than the QCD scale and 61 orders 
of magnitude smaller than the Planck scale. Obviously, all massive particles decouple the same way! 
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m\ = \/Ao ~ 10" 3 eV, is essentially given by the geometrical mean of the current value 
of the Hubble constant and the Planck mass, i.e. the smallest and largest energy scales 
conceivable in our universe: 

m A ~ V#o M P . (3.8) 

Eq. (^) provides a possible explanation for that. Moreover, Eq. (^) tells that the physics 
of the CC is naturally dominated by the set of sub-Planckian masses, irrespective of all 
the dynamical details of the low-energy fields with masses mi <C Mi, such as the SM fields. 
This idea completely frees the running of the CC from all kind of fine-tunings thanks to 
the smallness of our RG scale [i = H. In this suggestive scenario the running of the CC 
at any time is smooth enough, in particular also at the present time. At any epoch the 
rate of change of the CC is in the right ballpark to shift the value of CC in less than the 
value of the CC itself in that epoch. On the other hand, at higher and higher energies the 
RGE (|3.G| ) predicts a CC value increasingly larger. For instance, at the Fermi epoch, when 



— 1/2 

the temperature was of the order of the Fermi scale Mp = G F ~ 300 GeV, the Hubble 
parameter was of the order of H ~ T 2 /M P = Mp/M P and Eq. ([T|) predicts a typical 
value for the CC around A ~ H 2 Mp ~ Mp, which naturally fits with the value expected 
for the CC at the time of the electroweak phase transition. 

The origin of the Planckian mass operator on the r.h.s. of Eq. ( |3.6| ) could just be 
the indelible imprint left forever on the low-energy physics due to the decoupling of the 
sub-Planck mass fields just below the Planck mass scale. This permanent imprint may be 
thought of as a "relic" low-energy effect from the high energy dynamics of some fundamental 
RGE of the CC at the trans-Planckian scale \x > Mi, 

in which the functions F^fi/Mi) of the Planckian masses Mi depend on the underlying 
details of the trans-Planckian physics, e.g. string/M-theory. This Ansatz should hold good 
perhaps in the border line \i > Mp. Unfortunately we do not know the details of the RGE 
( |3.9| ) as we do not know the actual structure of the functions Fi(ii/Mi) for fi > Mj ~ Mp. 
Actually for fi 3> Mp this picture must break down as it probably does not even make 
sense to talk of the Hubble parameter because the metric need not to be the FLRW one. 
Indeed, for /i > Mp we just enter the realm of quantum gravity, where the metric itself 
is highly fluctuating. Hence the Ansatz ( |3.5| ) should be sensible only below the Planck 
scale. Then, and only then, we may set fi ~ H, and this identification should be better 
and better the smaller is the energy as compared to the Planck scale. It is only at these 
"low energies" that the soft decoupling of the sub-Planckian masses dominates the RGE. 
For instance, if the form factor takes the canonical form 

F(n/Mi) ~ F(H/Mi) = M2 H * H2 , (3-10) 

then for H 2 <C M 2 we may expand the term on the r.h.s. of (|3.9|) just to find 



Mf F(H/Mi) =M 2 H 2 -H i +^ + ... (3.11) 

i 
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In this way we arrive at some heuristic justification of Eq. (|3,l| ). At leading order in 
Eq. ( [3.9| ) reduces to our fundamental sub-Planckian operator in J3.6| ). 

At super-Planckian energies, H > Mi, the form factor ( |3.10 ) is of order one, and the 



r.h.s. of ( |3.9| ) behaves like Mf. In this regime we may expect an RGE of the form 

dln/j, (47r) 2 i 

Therefore, the typical CC at trans-Planckian energies, just in the upper neighborhood of 
Mp (// > Mp), becomes of the natural size M p and one may ask what to do with it. There 
is the attractive possibility that in this Planck neighborhood there is exact supersymmetry 
(SUSY), and if so there will be as many boson fields with mass Mj as fermion fields of the 
same mass, and since SUSY applies the sum on the r.h.s. of ( 3,12| ) could actually vanish 



identically. Then, at low energies we find that the CC is always controlled by our leading 
term H 2 Mp, and when there is a chance for the Mp contributions to appear, SUSY kills 
them automatically. So this would leave us with a well behaved CC at low energies, and 
all the dynamical details associated to the phase transitions below Mp (in particular the 
electroweak SM one) would be innocuous for the running of the CC . This property is 
robust within the low-energy regime (fi < Mp) and is guaranteed by the structure ( |3.6| ) 
of the RGE, independent of what particular speculation is made at fj, > Mp - e.g. Eq. 
Q3.12 ) and the aforementioned SUSY scenario. Of course the SUSY interpretation is only a 



possibility and we cannot be too conclusive. In fact, we cannot say much about the physics 
at trans-Planckian energies, not even at the border line /U > Mp, because the relation 



(|3.10[) is expected to be valid only for /i < Mp. However, this kind of situation is not too 
different from what we have with strong interactions in QCD. At high energies one meets 
asymptotic freedom, but in the infrared region the RGE of QCD tells us that the coupling 
constant grows. However one can not really conclude that it finally explodes because we 
are using an equation that is only valid in the perturbative regime. Similarly, here we 
ignore how the functions F([i/Mi) behave for [i > Mp and in particular the identification 
fi = H ceases to make sense, so strictly speaking we cannot use Eq. ( |3.10 ) to predict 



equation (3.12) for H > Mp. However, like in QCD, we can foresee a plausible trend 
in the behaviour for both the low and high energy regimes, specially if we invoke exact 
SUSY above Mp. Some more discussion on these issues, including the potential existence 
of non-local effects that might appear in the present approach, is provided in Appendix 1. 

4. FLRW cosmologies with a running cosmological constant 
4.1 Solving the model 

In the previous sections we have motivated our model. Let us now consider it in detail 
and show that it is useful and testable. One has to solve the coupled system of differential 
equations formed by Friedmann's equation, the equation of continuity (in the matter era, 
where pressure p = 0) and our RGE, i.e. the system formed by ( |2.3| ), ( |2.7| ) and ( |3.6| ) with 
H = H: 

f/A ' <rH 2 M 2 ■ .... (4.1) 



din H (4tt) 
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Here we have introduced the following mass parameter: 



M 



c,m;- 



(4.2) 



Furthermore, a = ±1 indicates the sign of the overall /?A-function, depending on whether 
the fermions (p = — 1) or bosons (a = +1) dominate at the highest energies. Notice that 
the mass Mi of each superheavy particle in (3.6) may be smaller than Mp and the equality, 
or even the effective value M > Mp, can be achieved due to the multiplicities of these 
particles. From Eq. fl3.7|) we see that the r.h.s. of (4.1) is of the order of the present value 
of the CC. 

Let us now eliminate the time variable and convert the equation of continuity 
into a redshift differential equation: 



dA dz da dp dz da 

dz da dt dz da dt ^ 







(4.3) 



Using the redshift definition ao/a = 1 + z and a = a H it immediately gives a very simple 
expression in which both a and H cancel out: 



dA dp 

dz dz 



3p 

1 + z 



(4.4) 



One can easily check that if A would not depend on the redshift (dA/dz = 0) then the 
previous equation integrates to pm a 3 = p° M Oq, i.e. we recover the old case. However, in 
general this is not so and now Eq. Q4.4j ) must be integrated together with (4.1) and (|2.3|). 
Eq. ( [4,ip can be transformed into a redshift differential equation by applying again the 
chain rule: 

. 2 dA , 2 dA dlnHdH 1 o dH 2 /A . 

(4vr) 2 — = (47r) 2 — — — — — — = -aM 2 — — . (4.5) 



dz 



d\uH dH dz 



dz 



To compute dH 2 /dz we recall Friedmann's equation ( |2 . 3| ) . Using the identity —k/a 
Hi n° K (l + z) 2 , it takes the form 



H 2 (z) 



8irG 



[p(z)+A(z)}+H 2 n° K (l + z 



(4.6) 



From this we have 



dH 2 8irG (dA dp\ rr2 ~ , 



Substituting Eq. (4^) into the previous equation and then the result into (|4 

3(1 -v)p(z) 



dp 

dz 



l + z 



K P ° C (1 + Z)=0. 



(4.7) 
we find: 
(4.8) 



Here p® is the critical density, and we have introduced for convenience two dimensionless 
parameters 

k = -2vVl\ (4.9) 
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and 

a M 2 , 

v = o . (4.10) 

12vr M 2 V ' 

Parameter k is related to curvature effects, and it is not independent of v once the spatial 
curvature k is known. Our model has one single independent parameter, v, which will 
play an essential role in the forthcoming discussions. The standard FLRW cosmology 
corresponds to v = 0. From now on v will parametrize all the cosmological functions 
that we obtain in our modified ( "renormalized" ) FLRW framework. For example, the one- 



parameter family of solutions of the differential equation(4 : 8) is completely analytical and 
reads as follows: 

p(z; v) = (p°m + p°) (i + z) 3{1 ~ v) - p° c (i + z) 2 • (4-11) 

The arbitrary constant has been determined by imposing the condition that at z = 
we must have p = p M . As we have said, the parameter k above introduces additional 
z/-effects due to non-vanishing spatial curvature. If we assume a > 0, then n > (resp. 
k < 0) corresponds to positively (resp. negatively) curved universes, i.e. closed (resp. 
open) cosmologies. For v = we also have k = and one recovers the expected result 
P = p\i (1 + z ) 3 ' i- e - po 3 = P%i a Q- However, for v ^ the parameter v really plays 
the role of a new cosmological "index" determining the deviations from the usual law of 
evolution with the redshift. Substituting ( [4.11| ) in (4.4) we may explicitly solve also for the 



^-dependent A, which becomes a function of the redshift: 

A(z; u)=A + p° M f(z) + p° c g(z) , (4.12) 

with 

f(z) = -^- Ul + zf 1 -")-!] , (4.13) 



Notice that the function f{z) is non- vanishing even if the spatial curvature is zero (k = 0), 
whereas g(z) ^ introduces curvature effects. To avoid confusion, we note that f(z) is 
well defined in the limit v — > 1. Similarly, the value v = 1/3 is non-singular in g(z). 

We have presented the CC and the matter density function as explicit functions of the 
redshift because it is the most useful way to present the result for astronomy applications. 
Eq. ( f4.1| ) can be trivially integrated with respect to H: 

A(z; u)=A + M 2 {H 2 (z; u) - H 2 ) =A + ^-M 2 (H 2 (z; v) - H 2 ) . (4.15) 

Computing dA/dz = (3 v / '8tt) Mp dH 2 / dz from this equation, and using ([4.7D and the 
equation of continuity (|4.4| ) to eliminate dA/dz, it is immediate to check that we are lead 
to Eq. ( [4.8|) . This shows the consistency of the whole procedure. Not only so; actually 
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Eq. ( 4.1 5| ) can be also useful from the astronomy point of view because it expresses a 
relationship between the CC and the Hubble parameter at any redshift. This correlation 
could be an experimental signature of this model, because it does not take place in the 
standard model. Furthermore, since we have already obtained the function A = A(z;v), 
we can use it in ( |4.15| ) to get the explicit function H(z; v). It reads 



h 2 {z]v) = Hi 1 1 + n 



.3(1-1/) 



1 - V 



+ 



i _ oo _ qo 
1-Zv 



(1 + zf - 1 - 2v 



3(1-*) 



1-V 



(4.16) 



It is easy to see that for v = we recover the standard result: H 2 (z; 0) = H 2 (z), where 



H\z) = H 2 {l + Qit [(l + *)3 _ i] + (i 

= Hi [n° M (i + zf + n° K (i + zf + n 



O 

A ) 



[(i + .) 2 -i]} 



(4.17) 



and therefore Eq. ( 4.16 ) constitutes a generalization of this formula for our model. The 
deviation should perhaps be testable (see below). 

Recall that the evolution in the remote past is obtained in the limit z — > +oo and the 



asymptotic evolution to the future corresponds to ; 



T. Then, some of the features 



of A(z; v) and p(z\ v) depending on the value and sign of the fundamental index u, are the 
following: 



For v < the CC becomes negative and arbitrarily large in the remote past, and 
at the same time the matter density infinite and positive, which is fine provided the 
latter dominates in the nucleosynthesis epoch. In the infinite future the CC becomes 
finite while the matter density goes to zero. For the flat universe the finite value of 
the CC in the asymptotic regime is positive and given by 



A(z 



-I) 



An + 



1 



Pm 



(4.18) 



This case is not incompatible with the measure of a positive CC in the recent past 
because all these models satisfy the boundary condition A(0) = An, and therefore if 
the CC is negative in the very early universe (anti-de Sitter space) it may just have 
changed sign recently. One can easily show that the transition redshift satisfies 



In 



ln(l + z) 



i + (i + fa) (n° /n°,) 



3(1+1/ 



(4.19) 



For example, for = 0.3, Q\ = 0.7 and v = —0.1 the transition from negative to 
positive A occurred around z = 1.7. This possibility cannot be excluded in the light 
of the present data, which are barely available at such redshifts. For smaller \u\, say 



5 This is strictly true only if the universe expands forever. If not, then it is approximately true, in the 
sense that it is valid near the turning point. 
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v = —0.05, the transition redshift becomes higher: z = 2.4. In principle, this v < 
case is not bad, except if \i>\ is too large, in which case the transition redshift would 
be too low and would have been detected. For example, for v = —0.4 the transition 
would be at z = 0.69 (see Section [| for more details). It is thus clear that the 
parameter v cannot be arbitrary and becomes restricted by experiment. Moreover, 
too large v (even if \v\ < 1) would also lead to problems with nucleosynthesis (see 
Section O). 



Another interesting case is to suppose that 

< v < 1 , (4.20) 

where the inequality signs are strict. Then the CC is infinite and positive in the 
remote past, and at the same time the matter density is also infinite and positive. 
Furthermore, the CC can be finite in the asymptotic regime while the matter density 
goes to zero, which is a double combination of welcome features. The value of the 
CC in the future is, in the flat case, 

A( z = -l)=A - T ^-p° M . (4.21) 

Of course this is just Ao if v = 0. However, for non-vanishing v the CC will be 
positive or negative in the asymptotic regime, depending on whether v < or 
v > Q\ respectively 6 . The change of sign from A > to A < in the course of the 
history of the universe can be of interest, see below. One can check that this feature 
could be maintained in the presence of the curvature term; in particular, this is so 
for k < if < v < 1/3, and for k > if 1/3 < v < 1. 

In the flat case, and for v > 1, the CC is finite and positive in the remote past, where 
it takes the value: 

A(z = +oo) = A + p° M . (4.22) 
However, this solution seems not to make much sense because the matter density 



(4.11) goes to zero in the remote past. This would not be the case for k > because 



then the density can go to infinity, due to the second term of (4.11). Notice that 



the flat case with zero matter and finite CC in the remote past could momentarily 
be considered as tenable in that perhaps the universe was first in a pure state of 
CC and then this CC decayed creating matter at much later times. Such "decaying 
CC" ( 4.22[ ) could be arbitrarily big if v — ► 1 + . Nonetheless, the v > 1 scenario (with 



or without curvature) has a big stumbling block: while the CC becomes large and 
negative in the infinite future, the matter density increases too. This is possible due 



to the balance of energy between A and p, Eq. (2.7), but a progressively more dense 
universe looks undesirable because does not seem to fit the trend of the observed 
evolution. 



6 In the A < case Eq. ( 1.21 ) can only be approximate because the universe will eventually stop expansion 



at some z in —1 < z < 0, see Eq. (2.8) 
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A critical case is v = 1. Here the matter density ( 4.11 ) becomes constant for flat space 
at all redshifts, and this does not look much sensible. This scenario, however, could 
perhaps be rescued for non-vanishing curvature as follows. For k > (corresponding 
to positive curvature k > in this case) Eq. ( [4.1 1| ) tells us that the matter density is 
infinite and positive in the remote past and very small in the long run future, which 
is what we want. Furthermore, this situation can be somewhat attractive because it 
smoothly matches up with case ( |4,20| ) for v — > l - . In this limit f(z) — > 31n(l + z) , 



g(z) -> {k/2)[z(z + l)/2 + 31n(l + z)] and Eq. fljlg) implies that the CC is positive 



and very big in the remote past, it eventually changes sign and it starts getting 
increasingly negative (remaining finite, though, see the previous footnote) in the 
future. It is easy to see that for sufficiently small (positive) curvature (which is in 
fact what we want in order not to depart too much from the flat case) , the transition 
into the negative CC regime will take place in the future at the redshift 

z ~ exp (-n° A /3n° M ) - 1 . (4.23) 

Thus for Q° M = 0.30, Q° A = 0.71 (implying Q° K = -0.01) we get z = -0.54, a far 
point in the future. We stress that although the v = 1 scenario is not possible in 
the strict flat space case, the choice of cosmological parameters that we have made is 
perfectly compatible with present day CMB measurements [Q, ||]. On the other hand 
the case of negative curvature would be a disaster because we encounter an infinitely 
negative mass density in the remote past. 

Some reflections are now in order. Take the flat case first. A most wanted situation for 
string/M-theory is, as we have mentioned in the introduction, to have negative (or zero) 
CC in the far future enabling to construct the asymptotic S-matrix states. As we have 
seen above, a necessary (though not sufficient) condition for this to happen in the present 
framework is to have v > 0. For instance, for v in the range ( [4.20| ) we can start with a 
large and positive cosmological constant in the early universe, which then decreases more 
and more and eventually, if v > J)^ , it becomes finite and negative in the asymptotic 
regime; in fact it can be rather large and negative if v — > 1~ . Actually it suffices that 
the CC is negative, no matter how small it is in absolute value, to secure the eventual 
stopping of the accelerated expansion and the disappearance of the event horizon. The 
largest possible value of v for which the CC can still change from positive to negative 
value is 7 v = 1. However, this limiting scenario is only tenable at the expense of having a 
positively curved universe. In contrast, the < ^ < 1 solution is in principle allowed in 
the flat case. We also remark that there is the possibility to have vanishingly small CC in 
the asymptotic future. This would be the case if v = f2^. The values of v for which A < 
in the asymptotic future are smaller than one, but since the present day estimate of the 
parameter is ~ 0.7 [||, the necessary v values imply a fairly large correction to some 
standard laws of conventional FLRW cosmology, especially in the flat case 8 . Whether we 

7 The v — » 1 + limit is troublesome because of the unwanted behaviour of the matter density in the 
asymptotic regime. Therefore, v = 1 sets a barrier and v > 1 should be considered unlikely. 

8 In the curved cases the solution cannot be obtained in closed form, and we shall not enter the details 
here. 
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can accept them or not is not obvious by now. However, if accepted, then it would hint at 
the "symmetry" approach to the old CC problem, in the sense that string/M-theory itself 
could perhaps provide that value of v as a built-in symmetry requirement. Nonetheless, 
before jumping to conclusions, we still have to check what values of v could be incompatible 
with nucleosynthesis. We do this in the next section. 

4.2 Restrictions from nucleosynthesis 

Needless to say, it is important to check what happens with nucleosynthesis in this model 
because a non-vanishing v may have an impact not only in the matter-dominated (MD) 
era, but also in the radiation-dominated (RD) epoch. We have seen that the index v enters 
the power of 1 + z in the expressions for p and A, and in the MD era we have (1 + z) 3 ^ 1- ^ 
rather than the standard behavior (1 + z) 3 . Similarly, in the radiation era one expects a 
behaviour of the form (1 + z) 4 ^~ u ^ . To check this we recall that in the RD era the equation 
of continuity ( |2.7| ) must include the p ^ term. For photons the radiation density pn is 
related to pressure through p = (1/3) pn, and we have 

A + p + AH p R = 0. (4.24) 

From the chain rule we can again trade the time variable by the redshift variable, and the 
previous formula becomes: 



This equation must now be solved in combination with ( |4.1| ) . We will not repeat the detailed 
steps. It is easy to see that the solution is obtained by simply replacing 3(1 — v) — > 4(1 — v) 
in Eq. ( 4.11| )-( |4.12 ). The radiation density at any redshift reads: 



P*M=[&+2=to) U + -) 4(1 ^ ) -^;U + ^ (4-26) 

where 

P ^2.5xl0- 5 / l0 - 2 / 9e (4-27) 

is the radiation density at present. In this case, since we are in a thermal bath of radiation, 
it is more natural to express the above result in terms of the temperature: 



Pr(T; v) = [Pr + 



n iM k /r^ 



Av ) \T J 2-Au \T f 



2 / rr \ 4f 



9*T 777 + 



30 " \T 2-Au 



j, \ 4(1-1/) f rjn \ 2 



(4.28) 



'-o ~ 2.75 K = 2.37 x 10 eV is the present CMB temperature. Of course for u — ► 
we recover the standard result 

PR(T)=p° R (^) 4 = ^5*T 4 , (4.29) 
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with = 2 for photons and = 3.36 if we take neutrinos into account. From these 
equations the restrictions imposed by nucleosynthesis are rather evident. Let us first of all 
quote the corresponding prediction for the CC in the radiation epoch, according to this 
model: 

A R (T; u) = A + p° M f R (T) + p° g R {T) , (4.30) 



with 



and 



f R (T) 



1 - v 



T 



4(1-") 



9R(T) 



T 2 



T 2 



2-4j/ 



T 2 



1 



4(1-*/) 



1 



(4.31) 



(4.32) 



Again we have separated the result into two functions fn(T) and Qr{T) in analogy with 
the MD epoch. We point out that the limit v — > 1/2 in equations ( 4.28| ) and (|4.3C| ) is well 
defined, as it was also the case with v = 1/3 in the matter epoch. 

In spite of the various possible scenarios that we have described before for the present 
MD era, we see that the RD era imposes additional conditions on the range of values of 
the cosmological index v: 



v = 1 becomes unfavored by nucleosynthesis. It is certainly ruled out in the flat case, 
otherwise the density of radiation at the nucleosynthesis would be the same as now 

1 scenario could still be argued if one 

Actually this 



- Cf. Eq. ( 14.281) . In the n > case the v 

accepts the law H oc T, instead of if oc T 2 , during nucleosynthesis, 
cannot be completely excluded on the basis of existing phenomenological analyses of 



the Friedmann equation in the nucleosynthesis epoch 51], but in general we shall 
stick here to the most conservative point of view. 

In the flat case, all u > 1 scenarios are troublesome because the density of radiation 
at the nucleosynthesis time falls below the one at the present time. For non- vanishing 
positive curvature, this situation could be somewhat remedied. However we already 
noticed that the v > 1 case was untenable because, irrespective of the value of the 
curvature, one predicts (in the MD epoch) a progressive growing of the matter density 
in the long run future. Therefore, the v > 1 case remains unfavored. 

On the basis of the most conservative set of hypotheses related to nucleosynthesis we 
conclude from the structure ( 4.28| ) of the radiation density, that the safest range for 
the v parameter is 

0<|i/|«l. (4.33) 

Both signs of v are in principle allowed provided the absolute value satisfies the 
previous constraint. In the following we will adhere to this possibility for most of our 
numerical analysis, although we shall leave open other possibilities for the theoretical 
discussion. 
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Furthermore, when comparing the relative size of the CC, Eq. ( pOl) , versus the 

at the time of the nucleosynthesis we naturally require 
For k = (flat space), the ratio between 



radiation density, Eq. ( 4.26 
that the former is smaller than the latter, 
the CC and the radiation density is fully determined by the index v 



PR(T) 



(4.34) 



Then it is clear that in order that the CC is, say, one order of magnitude smaller than 
the radiation density at the nucleosynthesis time, we must again impose Eq. ( 4.33; ). 



Then Eq. (|4.34|) leads to 



PR(T) 



f< 1. 



(4.35) 



It should be clear that the range ( |4.33j ) could already have been suggested from the 



behaviour of the matter density function ( 4.11 ) alone, if we are not ready to tolerate 



a departure from the exact (1 + z) law in our MD era. However, let us stress that 
in our era there is no crucial test (at least an obvious one) emerging from the present 
values of the cosmological parameters that is sensitive to the deviations v ^ in a way 
comparable to the nucleosynthesis test. Maybe the high precision future experiments 
(see Sections [| and ^) can put a remedy to this. At the moment the restriction on 
v coming from nucleosynthesis alone coincides with our general will to remain in the 
framework of the effective field theory approach introduced in Section |3|. Indeed, Eq. 
( |4. 33| ) implies - see Eq. (4.10)- that the effective mass scale M cannot be much 
larger than Mp, and in particular M < Mp is a natural possibility. This is the 
kind of picture that we wanted from the general discussion of our RG framework in 
Section ||. So, indeed, the restriction from nucleosynthesis tells us that we cannot 
play arbitrarily with the value of the new cosmological index v, if we want to get a 
consistent picture both theoretically and experimentally. 



The essential issue is whether the restriction ( 4.33 ) leaves still some room for useful 
phenomenological considerations at the present matter epoch. The answer is that it does. 
Assume for definiteness that a = +1, and let us take the most natural value for the 
cosmological index v: 

2.6 x 10~ 2 , (4.36) 



1 



12 tt 



which corresponds to M = Mp in ( 4.10 ). Obviously it lies in the nucleosynthesis suggested 
range ( |4.33| ). Let us next circumscribe the calculation to the flat case, where g(z) = 0. By 
expanding the function f(z), Eq. ( 4.13] ), in powers of v we immediately find, in first order, 



A(z)~Ao + up° M [(l + z) 3 -l] . 
What about the numerical effect? It is given by 



An 



(4.37) 



(4.38) 
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Taking n° A 



0.3, n° A 



0.7 and z = 1.5 (reachable by SNAP ||) we find <5 A = 16. 



7o, 



namely a sizeable effect that should be perfectly measurable by SNAP. For values of order 
v = 0.1 the previous correction would be as large as 47.8%. This v value would correspond 
to an effective mass scale (4.2) of M < 2Mp. 



4.3 The ^-dependent Hubble parameter 

It is also useful to realize that the value of the Hubble parameter at a given redshift, is 
different in our model with respect to the standard model, see Eq. ( |4.16| ). To check whether 
this deviation is testable in the near future, let us see how much Eq. (|4.16| ) departs from 
Eq. ( 4.17 ). Since j/ < lwe can expand Eq. ( 4.16| ) in powers of this parameter and subtract 
the standard result ( 4.17 ). For simplicity let us take the flat case (n = 0). The relevant 
subtraction is, in first order of v, 

A(z;u) = H 2 (z;v)- H 2 {z) ~ - v H$Q,° M {l + (1 + zf [31n(l + z) - 1]} . (4.39) 

Therefore, the relative deviation of the Hubble parameter in our model with respect to the 
standard one is, at any given redshift z, the following: 



5H(z;v) 



1 A(z;u) 
2 



H(z;u) - H(z) 

W) 

_1 o0 l + (l + z) 3 (31n(l + z)-l 

2 uUm i + n% [(i + zr-i] 



(4.40) 



where H(z) is the standard value given by the square root of ( |4.17| ). We see that the 
correction is negative for v > 0, in which case the model predicts a Hubble constant smaller 
than expected at any redshift - or larger if v < 0. Of course we have 5H(z = 0; v) = for 
any u, as expected, because at the present time we have input the values of the cosmological 
parameters. Let us check numbers for the v > case by considering the future SNAP 
experiment (see Section || for more details). Take a very distant supernova at z = 1.7, 
and assume that v is given by Eq. (|4.36| ), and that £1° M = 0.3. Then Eq. ( [4.40 ) gives 
5H ~ —2.4%. The deviation is not big. Fortunately, the other cosmological parameters 
are much more sensitive, and our simulation analyses in Section |?j will not be dependent 
on H. For both signs of v the Hubble parameter H(z; v) tends to a constant determined 
by the cosmological term in the infinite future (z — ► —1), as in the standard case. So we 
are just testing the different rate at which H (z; v) goes to that constant as compared to 
H(z;0). Let us find the asymptotic future regimes for small v, even in the presence of 
curvature. From Eq. ( 4.16| ) the result is 



H\z 



-l;i/) 



m 



K 



1 



^H 2 [ni-u(n° M + n° K )] 



(4.41) 



and for the standard model is of course 



H 2 (z = -1;0) 



(4.42) 
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Therefore, in the asymptotic regime the relative difference is, in first order in v, 

1 oo , no 
8H(z = -1; v) ~ -- 1/ JU ^ . 



(4.43) 



In the flat case (fi^ = 0), with as in ( 4.36 ) and the usual values for the cosmological 
parameters, it gives 5H(z = — 1;z/q) ~ —0.5%. 



4.4 The ^-dependent cosmological sum rule 

From the Friedmann Eq. ( f4.6| ) we immediately find, for any v, 

£Im(z; v) + il\(z; v) + tt K (z; v) = 1 , 

where 



^m(^) 



Pm{z;v) 
Pc{z;v) 



Q.\(z; v) 



Here 



A(z;i/) 

p c (z; ^) 



-k 



H 2 (z;u)a 2 



Pc{z;v) = 



! 7T Gtv 



(4.44) 



(4.45) 



(4.46) 



is the value of the critical density at redshift z. Eq. ( 4.45j ) reduces to the standard one for 
v = 0. The sum rule ( 4.44 ) is exact at any redshift and any value of v because Hubble 's 
"constant" and the critical density correspond also to that redshift. Now, having said 
that, consider the flat case which is well motivated by inflation. The sum rule is then 
Qm(z;u) + Q.^(z;v) = 1, and can (in principle) be checked experimentally by measuring 
the Hubble constant, the matter density and the CC at a given redshift and then computing 
the values of the first two parameters in ( 4.45| ). However, suppose that we did not suspect 
that the CC itself is a function of the redshift and we just assumed that A = Ao for all z. 
Then it is clear that the previous sum rule would fail, because the evolution of the CC and 
matter density are indeed related. 

Let us consider the explicit expressions for the cosmological parameters in the flat case. 
Using the previous formulae for CC, matter density and Hubble parameter in our model, 
we find 



8irGp M {z; v) 
3H 2 (z;u) 



ni r (i + z ) 3{1 - u) 



1 _l_ M 
1 "i" 1-v 



(1+*) 



3(1-0 



(4.47) 



and 



£Ia(z;v) 



8ttGA(z;v) = SA + lh n°/ [(1 + *) 
3H 2 (z;u) 



3(1-0 



1 ^ \-v 



[1 + z 



,3(1-0 



The flat space sum rule is clearly borne out at any z and for all v, 

U M (z; v) + Q\(z; v) = 1 , 



(4.48) 



(4.49) 
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after making use of the present day sum rule. It is also easy to check the following be- 
haviours in the remote past and in the infinite future: 



Qm(z = oo) = 1 — v , Q\(z = oo) = V , 

n M (z = -l) = 0, n A (z = -l) = l. (4.50) 

Clearly the flat space sum rule is satisfied also in both of these extreme (past and future) 
regimes: 

Q M (z = oo) + Q A (z = 00) = I = Vt M {z = -1) + tt\(z = -1) . (4.51) 

In particular it is interesting to note that in the infinite future z — > —1 we have Q A (z; v) — ► 1 
and Qm(z] u ) ~ * both in the standard model (u = 0) and in our model, irrespective of 
the value of v. Indeed, the point (£Im = 0, Vt\ = 1) is a fixed point to which the cosmic 
flow is attracted to, and this occurs even for nonzero v. However, we remark a difference 
in behaviour in the remote past: in the standard case Sl A (z — » 00) — > whereas in our 
model it never vanishes and it tends asymptotically to £l A (z — » 00) — > v. In the usual case 
at hand, with v given by Eq. ( [4.36 ) and a flat scenario, we find that the $7 a parameter in 



the remote past was £l A (z — ► 00) = 0.026, i.e. some 27 times smaller than it is now, and it 
never vanished, contrary to the standard FLRW model. 

With some more algebra one can also check the fulfillment of the ^-dependent cosmo- 
logical sum rule in the curved case, for any n. Let us show it for z = —1 (the asymptotic 
regime), as it can be illustrative. From Eq. ( pl| ) it is clear that p{-l;v) = even for 
k 7^ 0, so fi(— l;v) = too. On the other hand from the complete expression ( 4.12j ) at 
z = —1 (and no expansion in v) we have 

A(-l; u)=A - I ^-p° M - T ^- (p° c - p° M - A ) . (4.52) 
Moreover, from the ^-dependent Hubble's constant ( |4.16| ) we have 



1 



O 4- O 
1-1/ 



(4.53) 



where one could also have used the first line of ( 4.41 ). Finally we bring together these two 
pieces and get 

n / ^ 8ttGA(z = -1) 

n A { z = -1 = — = 1, 4.54 

3H z (z = —1; v) 

so the test is successful for all v and arbitrary curvature, in the case z = — 1. For arbitrary 
z, k and is, the sum rule also holds but it is more cumbersome to check. 

4.5 Predicted deviation of $7a(z; v) from the standard Sl\{z) 

It should be useful to compute the deviation of the z- and ^-dependent £l\ parameter in 
our framework and in the standard FLRW case. Let us define the deviation parameter: 

5n A {z; v) = , (4.55) 

U A (z;0) 
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where Q\(z;0) = ^a(^) is the standard one at redshift z. It is clear that the deviation 
( 4.55[) must increase more and more with redshift because the independent parameter does 
not go to zero in the past, in contrast to the standard £l\(z), as shown by equation ( 4,5C| ). 
In the following we restrict the evaluation of ( f4,55| ) to the flat case. Substituting Eq. ( 4.4^ ) 
in (4.55) and performing an expansion of the result in powers of v we find, in first order, 



+ 



1 + 30° 



M 



[l + z) 3 ln(l + z) 



Notice that this formula vanishes exactly at z = 0, 



(4.56) 



50 A (0; v) = 0. 



(4.57) 



as it should because our model and the standard one with the z-independent CC have the 
same initial conditions at z = 0. Moreover, Eq. (4.56) has the two expected limits for the 
remote past and future: 



5S1a(oo; v) 



oo . 



Sn A (-l;u) = 0. 



(4.58) 



In order to visualize the numerical impact of ( 4.56| ) in high-z supernovae tests, let us use 



the usual set of inputs: v as in flOp , n° M = 0.3, Q° A = 0.7, and take z = 1.5. The 



result is JJ1a(1.5; u) = 20.4%. Certainly a 20% correction should be measurable at SNAP. 
Already for z = 1 the result is sizeable: 6£L\(l;v) = 10.2%. This 10% effect could be 
considered as the relative deviation undergone by Q\ when we compare the known set of 
high redshift supernovae used in this paper, whose average redshift is around z = 0.5 (and 
from which the parameters O^ = 0.3, Q A = 0.7 were determined), with the central value 
z = 1.5 of the highest redshift set to be used by SNAP. 



4.6 Deceleration parameter 

It is also interesting to look at the deviations of the deceleration parameter with respect 
to the standard model. It is well-known that there are already some data on Type la 
supernovae located very near the critical redshift z* where the universe changed from 
deceleration to acceleration [^]. But of course the precise location of z* depends on the 
FLRW model and variations thereof. In our case z* should depend on our cosmological 
index u, i.e. z* = z*(y). The definition of deceleration parameter leads to 



a 



aa 

{Z]U) = ~~dJ = ~ aR\z-v) ~ 2 



- [Q,m{z\v) -2Q,a.{z;v)\ . 



Equivalently, 



(4.59) 



. . H 1, . 1 dH 2 (z:u) 

,(,;,) = -l-- 2 = -! + -(!+,) 



(4.60) 



For simplicity in the presentation, let us consider the flat case. Substituting either Eq. 
(4.47) and ( |4.48j ) in ( [4,59| ), or just Eq. (4.16) in (4.6C), and expanding in first order of v 
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we find the ^-dependent deceleration parameter, 



q(z\v) 



1 



2 1 + n° M [(i + 



(4.61) 



31n(l + z) + 



n° M [(1 + z) 3 - 1 - 3(1 + z) 3 ln(l + z) 



i + n° M [(i + 



i] 



For z 



we recover the standard FLRW result 

3 



1(0; u) 



-1 + 2°°M 



0° 



o A 



?0 



(4.62) 



for a flat universe. Of course the last formula is independent of u, because we have 
normalized our inputs to reproduce the cosmological parameters at present. On the other 
hand for v = 0, but any z, 

3 n° M (i + z) 3 



l(z;0) 



-1 + 



2 1 + 0° [(1 + 



1]' 



(4.63) 



This is the deceleration parameter as a function of the redshift for a standard flat FLRW 
universe. This result vanishes at the redshift z = z* , where 



-1 + 



09 



(4.64) 



for n° A 



0.7 and £l° M 



0.3. Hence the value ( [4.64 ) represents the transition point 
from a decelerated regime q(z\v) > (corresponding to z > z*), into an accelerated 
one q(z;v) < (corresponding to z < z*), within the flat FLRW standard model. It 
should be clear that this transition, for any curvature, does not represent the border 
crossing from a matter dominated universe, pm{z) > A(z), into a CC dominated universe, 
Pm(z) < A(z) (actually this crossing occurs later at z ~ 0.33); it rather represents (see Eq. 
Q4.59D ) the transition from the era where £Im(z) > 2 0a ( z ) (decelerated expansion) into 
the era where VLm{z) < 2fl\(z) (accelerated expansion) - therefore from pm{z) > 2A(z) 
to pm(z) < 2A(z). For v ^ 0, and in the particular case of the flat space, it defines the 
transition from Q\(z;v) < 1/3 to £l^(z]v) > 1/3. Then it is not difficult to see that the 
following inequality defines the value of z*(y): 



1 



(1 + *) 



O 
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< V 



:i + z) 



l + -ln(l + z) 



1 



(4.65) 



If this inequality is satisfied, it means acceleration; if it is violated it means deceleration. 
The ^-dependent transition point z*(v) is defined by the equality of both sides. Notice 
that for v = it immediately reproduces the previous result (4.64). The inequality cannot 
be solved analytically for v ^ 0, but in the next section we provide the numerical results. 
Already analytically it is obvious that for v < the critical redshift z*(v) will become 



smaller (closer to our time) than ( 4.64 ), whereas for v > the value of z* will be larger, i.e. 
the transition from deceleration to acceleration occurs earlier. In the presence of curvature 
(k, ^ 0) the analytical expressions defining the transition point are more cumbersome and 
we limit ourselves to present the numerical results in Section B 
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4.7 A note on inflation and the RG approach 

As a special theoretical issue concerning the RG framework presented here, let us say a few 
words on how to potentially incorporate inflation. Our ^-dependent cosmological equations 
(Cf. Sections 4.1 and 4.2) predict, for < v <C 1, that at higher and higher energies there 



is a simultaneous increase of both the matter /radiation energy density and CC term. For 
example, at some Grand Unified Theory (GUT) scale M x , where H ~ M^/M P , our RGE 
(|4.l| ) naturally predicts A ~ M^- During the RD epoch the CC is always smaller than 



the radiation energy density by a factor u <C 1, Eq. ( 4.35 ). However, this picture must 
break down at the very early epoch where radiation is not yet present (matter-radiation 
are still to be "created"). It is conceivable that this fast inflation period occurs near the 
Planck scale, following e.g. the anomaly-induced mechanism suggested in [fiS|| , which is 
a modification of the original Starobinsky's model pSf] . If there is SUSY, as speculated 
in Section |3|, and this symmetry breaks down at some energy near Mp, then there is no 
contribution to the CC running above that energy. However, just below Mp a CC of order 
M p is induced (Cf. Eq. fl3.12| )) due to the mismatch between the boson and fermion masses 
at that scale, and so inflation can proceed very fast. As inflation evolves exponentially the 
scale fj, = H decreases and the SUSY particles decouple progressively. Since the total 
number of scalar and fermion d.o.f. lessens with respect to the number of vector boson 
d.o.f., the anomaly- induced inflation mechanism becomes unstable and it finally leads to a 
FLRW phase- see the details in Ref . |49| . At this point the RD epoch of the FLRW universe 
starts: the radiation is supposed to have emerged from the decaying of that vacuum energy 
density, of order M p . Of course we have H < Mp after inflation, and the RGE has already 



changed to Eq. (|4.1|). So, following the above discussion, we are left again with the ratio 



(4.35), which insures <C pp and hence safe nucleosynthesis. The details of the combined 
mechanisms will not be discussed here, but it is clear that the model of Ref. |4{| can be 
naturally invoked in our CC approach because that model is based on the decoupling of 
the heavy degrees of freedom according to the RG scale ([T^), exactly as in the present 
framework. 

5. Numerical analysis of the model 

In order to see the behavior of the most representative parameters describing the universe, 
we analyze numerically the results obtained in Section || for the physically interesting values 
of the cosmological index, \v\ <C 1. In particular we use v = 0, ±z^o, ±2vo, ±0.1, where uq 
is given by ( 4,36) . 



We will first of all concentrate on the flat case and later on we consider an extension 
to k = ±1 . In the following we take £1® M = 0.3 and = 0.7 at z = for a flat universe, 
U° M = 0.3, tt° A = 0.5 for an open universe, and tt° M = 0.4, tt° A = 0.8 for the closed 
case. Spatially curved universes are not favored nowadays by CMB data [Q], but we would 
like, nevertheless, to show their behaviour for some cases which significantly deviate from 
flatness. The kind of study we present here is mainly based on supernova data, and we treat 
this analysis independently from CMB measurements. Implications of the RG framework 
for the CMB will not be discussed here. 
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5.1 Flat universe 



Let us start with an universe with flat spatial section. In this case the evolution of the 
matter density and of the CC is shown in Fig.|l]a,b. These graphics illustrate Eq. ( 4,11| ) 
and ( |4.12| ) for k = 0. As a result of allowing a non-vanishing /^-function for the CC 
(equivalently, v ^ 0) there is a simultaneous, correlated variation of the CC and of the 
matter density. The evolution of pm(z; v) and A(z; v) with z and v is very relevant because 
these functions appear directly in the luminosity distance expression - Cf. Section [6|. 

Comparing with the standard model case v = (see the Fig.[j]a,b), we see that for a 
negative cosmological index v the matter density grows faster towards the past {z — > oo) 
while for a positive value of v the growing is slower than the usual (1 + z) 3 . Looking 
towards the future {z — > —1), the distinction is not appreciable because for all v the 
matter density goes to zero. The opposite result is found for the CC, since then it is for 
positive v that A(z; v) grows in the past, whereas in the future it has a different behaviour, 
tending to different (finite) values in the cases v < and < v < 1, while it becomes 
— oo for v > 1 (not shown). We have made some general comments on these behaviours in 
Section |3|, and given the limiting formulas for these cases; here we just display some exact 
numerical evolutions with z within the relevant intervals. 




Figure 1: (a),(b) Future and past evolution of the matter density pm{z;v) and the cosmological 
constant A(z; v) for a flat universe (k — 0) and for different values of the fundamental parameter 



v of our model (v$ is defined in Eq.(4.36)). In both cases fl° M — 0.3 and £1^ = 0.7. The solid line 
represents the standard model case (A = const.), whereas the various kinds of dots/dashes represent 
different amounts of evolution of A. 



In the phenomenologically most interesting case \v\ < 1 (see Section L2) we always 
have a null density of matter and a finite (positive) CC in the long term future, while for 
the far past yields A = ±oo depending on the sign of v. In all these situations the matter 
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density safely tends to +00. One may worry whether having infinitely large CC and matter 
density in the past may pose a problem to structure formation. From Fig.[l|a,b it is clear 
that there should not be a problem at all since in our model the CC remains always smaller 
than the matter density in the far past, and in the radiation epoch z > 1000 we reach the 
safe limit ( |4.35 ). Actually the time where A(z;v) and Pm( z ', u ) become similar is very 
recent. Take, for example, the flat case and assume the usual values of the cosmological 
parameters as in Fig. HI: then for v = (—2vq, —vq, , +v$, +2vq) equality of CC and matter 
density takes place at z = (0.29,0.31,0.32,0.34,0.36) respectively. For larger values of u 
(still in the \v\ <C 1 range), say v = (—0.1, +0.1), we find z = (0.27,0.43). In all cases the 
equality of matter density and CC corresponds to very recent times, and so the evolution 
of the CC in this model never prevented structure formation. 




Figure 2: £Im{z\v) and £l\(z]v) versus redshift for a flat universe and different values of the 
parameter i/, with VL° M and as in Fig.KJ. 

Related to the evolution of pm{z\v) and A(z; v) are the cosmological parameters 
Q,m{z;v) and £l^{z]v) respectively - see Fig.^a,b. They are sensitive to the evolution 
of the corresponding energy densities and at the same time to that of the Hubble expan- 



sion rate H(z;v) - see Eq. ( 4.16 ). The evolution of £Im{z;v) and Q,\{z]v) tells us how 



the present values and Q\ differ from the corresponding values in the past and in the 
future for the standard model case {y = 0) and the present model case [y 7^ 0). Although 
Q\ is found to be at present of the order of 70% of the matter/energy in the universe, as we 
approach z ~ 1 its contribution was only a quarter of the total. The exact value obviously 
depends on the value of v (Fig.|2|b), being already ~ 10% larger /smaller than for a con- 
stant A for v = ±z^o at redshift z = 1. If we go further back in time, fl\(z — » +00) always 
diminishes, and tends asymptotically to while in the standard case £l\(z — > 00) = 0. 
Just the opposite occurs for Qm(z; v) (Fig.|2|a), since we are showing the flat case and so 
the sum has to be 1. 
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Figure 3: (a),(b),(c) Relative deviations of A(z;i/), f^z;^) and H(z;v) from the standard CC 
model (f = 0) versus redshift for different values of v; (d) Rate of change of A(z; v) (normalized to 
the present critical density) with redshift and for various v' s. In all cases we consider a flat universe 
with the same parameters as in the previous figures. 



In Fig.|3]a,b,c we show the three deviations of the parameters <5A, b^l^ and 5H with re- 



spect to the standard model case (y = 0), as defined in Section 4.2, 4.5 and [43| respectively. 
Consider first the deviation of A(z; v) caused by the running. As a function of z and u, 
A(z; is) is very sensitive to v at large z. Thus at z ~ 2 the increment of A is of 50% 
for v = 2fo, and once again the effect is higher for negative v (Fig. |3]a). Similar or even 
larger numbers are obtained for 5$lf\_, which attains e.g. 70% under the same conditions. 
In contrast, the deviations 5H of the Hubble parameter from the standard value are much 
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Figure 4: Evolution of the deceleration parameter, q, for a flat universe (k = 0) with f2j[ f = 0.3 
and = 0.7. We see that the transition from a decelerated (q > 0) to an accelerated (q < 0) 
universe takes place earlier in time (larger redshifts) for larger v > (see text for the exact values) . 




Figure 5: Evolution of the matter density for an open universe with = 0.3 and £1° = 0.5 (left) 
and for a closed universe with fi^ = 0.4 and = 0.8 (right). Differences in the curves are more 
dependent on the value of fl° M than on the curvature. 

smaller (Fig.||c), around 6%. Finally, the CC variation rate with redshift (normalized to 
the current critical density), (l/p^) dA/dz, is presented in Fig.||d. These curves are non- 
symmetric in the sign of v when z becomes large. The effect is more important for negative 
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v. As one can already see from A(z; v), Eq. ( 4.12| ), as we go further in redshift the running 



increases more quickly for v < than for v > 0. 

To conclude the analysis of the flat case we consider another relevant exponent describ- 
ing how the universe evolves, the deceleration parameter q. This one is fully sensitive to the 
kind of high-z SNe la data under consideration. The ^-dependence of this parameter was 



discussed in detail in Section 4.6. The transition point between accelerated and decelerated 
expansion is a function of v : the more negative is v, the more delayed is the transition 
(closer to our time)- see Fig.|I|. If v > 0, the transition occurs earlier (i.e. at larger z). 
While in the standard case, and for a flat universe, the transition takes place at redshift 



0.67 -Eq. ( 4.64 ) -, it would have occurred at z = 0.72 and z = 0.78 for v = vq and v = 2vq 
respectively, and at z = 0.63 and z = 0.59 for v = and v = —2vq (Cf. Fig. |j). For 
v = (—0.1, +0.1) the effect is quite large, namely the transition would be at z = (0.53, 0.91) 
and hence there is a correction of (—21%, +36%) with respect to the standard case. 

5.2 Curved universe 

For the small values of v that we present in the previous figures, the differences between 
a flat k = universe and a k = ±1 one are not evident. A trivial variation is that 
coming from the different choice of the present-day values of and f2 A . Besides, we 
have more marked ^-dependent features than in a flat universe, but the main character- 
istics remain the same. In order to see the differences for the curved case, we show in 
Figures ||| and || the most representative parameters: pm, A, f^A and dA/dz, both for 
positive and negative curvature. Qualitatively, the behaviours are similar to the flat case. 
However we note that it is for positive curvature, i.e. closed universe, that we find the 
most dramatic numerical differences with respect to the flat case for each value of the 
cosmological index v. In particular, for closed universes we observe a faster running of 
the CC, especially for v < (this fact is confirmed when adopting other sets of cosmo- 
logical parameters for closed universes). Therefore, k = +1 universes would represent 
the most favored case for the possibility of the observational detection of the CC run- 
ning in our model. On the other hand the open universes differ numerically very slightly 
from the flat case even though the degree of (positive or negative) curvature chosen in 
the two examples shown in Figures [| and ^ is the same, namely Qk = ±0.2. As already 
advertised, from the point of view of the CMB data the curved cases under considera- 
tion would be excluded because Qo = 1.02 ± 0.02 |5j. However, we should still be open 
minded to the possibility of non-flat universes and maintain full independence of the two 
sets of data, CMB and high redshift SN la. Although we do not display the behavior 
of the quantities presented in the previous section due to their similarities, it might be 
interesting to comment where the transition from deceleration to acceleration takes place. 
In the case of the open and closed universes defined above the transition redshifts for 
v = (-0.1-^o,-2^o,0,i/o,2^o,0.1) read: z = (0.36,0.42,0.45,0.49,0.55,0.61,0.77) for 
open, and z = (0.49,0.53,0.56,0.59,0.62,0.65,0.73) for closed. So, the width of the red- 
shift interval is almost the same as for the flat universe, but the transitions tend to occur 
at lower (higher) redshifts for open (closed) universes. 
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Figure 6: Cosmological constant (top), density of cosmological constant (middle) and derivative 
(bottom) evolutions for an open (k < 0) universe with Q° M = 0.3 and = 0.5 (left) and for a 
closed (k > 0) universe with = 0-4 and = 0.8 (right). Numerical changes with respect to 
the flat case are see to be most pronounced for the closed case. 
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6. Magnitude-redshift relation and Type la Supernovae 



The analysis of supernova magnitudes allows to test different cosmological models since 
magnitude depends on the dynamical evolution of the universe. This is true not only for 
supernovae but for all standard candles, i.e., objects whose absolute magnitude M (or 
intrinsic luminosity C) we know, and whose apparent magnitude m (or received flux J 7 ) 
can be measured at a given redshift. For cosmological purposes related to the expansion 
history of the universe, SNe la are the best candles since their luminosity enables detection 
up to very high redshift. Besides, although there are no standard candles in nature, in the 
SNe la case the variety is well accounted for by a tight correlation between magnitude at 
maximum and decline of the light curve [Q, |36|, and the candle can be calibrated through 
such correlation. A way to parameterize this effect is through the stretch factor, s, used 
by the Supernova Cosmology Project (SCP) [§]. The stretch factor method expands or 
contracts by a factor s the time axis of every supernova light curve to fit a fiducial one. The 
stretch-corrected SNe la magnitudes of the sample are then fitted to obtain the cosmological 
parameters. 

The apparent magnitude obtained at different redshifts is related to a given cosmolog- 
ical model via the magnitude-redshift relation. One starts from the notion of luminosity 
distance, di, related to the received flux T and the absolute (intrinsic) luminosity L through 
the geometric definition Q: 

F = ^F L - ("J 

Then the logarithmic relation between flux and (theoretical) apparent magnitude reads 



rn 



lh 'z,H ,n M ,n A ) = m + 5i ogl0 [H d L (z,H ,n° M ,n° A )} . (6.2) 



In the last equation, terms have been defined in order to collect all the dependence on the 
current value of the Hubble parameter Hq into the expression 

M = M - 5 log 10 H + 25 . (6.3) 

This way all the model dependence is encoded in the luminosity-distance function di = 
d^iz, Hq, f2^, 0, A ). Notice that the combined expression Hq di entering the argument of the 
log on the r.h.s. of Eq. ( |6.2| ) is Hubble constant-free. On the other hand, from the FLRW 
metric ( |2.6|) , the luminosity distance of a source at (dimensionless) radial coordinate r and 
redshift z is given by di = oq (1 + z) r, where ao = a(*o)- So we need to compute r as a 
function of the cosmological parameters. Since our FLRW universes have been modified by 
the renormalization group effects represented by the z^-parameter, the luminosity-distance 
relation takes a slightly different form as compared to the standard one [JIJ. This can 
be foreseen from the generalized structure of the ^-dependent expansion rate ( [4.16 ) or, 



equivalently, the ^-dependent cosmological constant parameters ( 4.45| ). It means that in our 



modified FLRW model the luminosity distance becomes a function of z parameterized by v 
and the present day values of the cosmological parameters: d^ = di(z, Hq, O^-, Q a ; v). For 
v = this function reproduces the standard result. The explicit derivation of this function 
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for v ^ follows steps similar to the conventional case, namely one starts considering the 
equation adr = dt a/1 — k r 2 for a null geodesic along a radial direction, which follows from 
the FLRW metric ( |2.6[ ). This can be rewritten as 

dr 1 + z , 1 dz ,„ „. 

= dt = — - , (6.4) 

yl-kr 2 a a H(z) 

where H{z) is the expansion rate at redshift z. Upon integration on both sides we have 

dr' 1 l' z dz' 

1 (6.5) 



o v 7 ! - kr' 2 «o Jo H(z',n° M ,n° A ; 



where we have taken into account that in our case the expansion rate O^, fi9 ; z/) is 
the ^-dependent function H(z; v) given by Eq. ( 4.16| ). After trivial integration of the l.h.s. 



of Eq. (|6.5| ) for k = 0, ±1 one finds the desired radial function r = r(z,k, Hq,Q° m f2°). 
Trading the curvature parameter k for Vt Q K = —k/H^ a^, one immediately finds the exact 
luminosity-distance function 



with 



sinx, VP K < 
^(x) = I x, n° K = (6.7) 

sinhx, > . 

Here the difference with respect to the constant CC case is encoded in H (z, Q° M , £1%; v). 
For v = 0, H(z, ft° M , $7° ! ^ = 0) becomes the standard FLRW function ( [4.17|) , and the 
luminosity distance (|6.6|) also reduces to the standard form. 



We use the magnitude data from the SCP (Supernova Cosmology Project) ||. The 
set includes 16 low-redshift supernovae from the Calan/Tololo survey and 38 high-redshift 
supernovae (Figj7|) used in the main fit of 0. Observational data are corrected such that 
they can be used as the intrinsic magnitude of the object. The effective value of the 
magnitude is obtained from that at the peak of the light curve according to the expression: 

m e Jf = m p ™ k + a{s - 1) - K BX - A x = m eff + a{s - 1) . (6.8) 

Here a is the parameter describing the correlation between maximum brightness and rate 
of decline of SNe la; s is the stretch factor mentioned above; X is the observed band; Kbx 
is the correction to the change from the emitted B-band to the received X-band, and Ax 
is the galactic extinction. 

Figure [7] represents, in the magnitude-redshift space, the data obtained by the SCP 
Collaboration together with the predicted magnitude-redshift relation for our model based 



on the scale (3J5) and two RG models based on the scale fl3.3|) [23, 29]. At low redshift 
all the models are equivalent as it occurs with all the alternative models to the standard 
one. At high redshift they display small variations with respect to the "constant CC" 
cosmological model. As it is seen graphically, current data are not able to distinguish 
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Figure 7: The magnitude-redshift relation is similar for small variations of the cosmological con- 
stant in different RG models in the literature, e.g. models with very light neutrinos and scalar fiels 
that could change the sign of the [3\ function |2^, , or the opposite case involving the heavy fields 
of the SM (^9). Also shown is the standard CC case and the present model for v = ±2vq. Points 
represent the 54 SNe used in the main fit from Q 

between models. At present, this kind of models can only be favored theoretically (see 
Section ^ for the future prospects). 

In order to determine the cosmological parameters we use a x 2 -statistic test, where x 2 
is denned by the difference between the theoretical apparent magnitude and the observed 
one: 

{m + 5 log 10 (H d L ( Zi ,H , M , Oy, v)) - mf f - a{ Si - l)} 2 
X 2 (0 M , Ov, v, a,M) = \ ± 2 

i 

(6.9) 

As we are only interested in the cosmological parameters associated to the acceleration of 
the universe, we marginalize over a and M |2|, [541. That means that we minimize x 2 by 
integrating over all possible values of a and A4. 



6.1 Results from current data 

In this section, confidence levels surrounding best fits (minimum \ 2 ) are given by contour 
lines of constant x 2 > which represent la, 2a and 3a confidence levels respectively. 

A distinguishing feature of all the models is the evolution of the cosmological constant, 
thus we can make a first test in order to see whether such an evolution is consistent with 
current data. We adopt a generic form for the Hubble parameter which describes any CC 
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running: 



H\z) 



Hi 



O 

M 



(1+*) 



+ 



1 dA 



p: 



,0 



r/2 



z + (1 + z) 



(6.10) 



This equation parameterizes generically the deviations from the standard law ( [4,17| ) through 
a series expansion of A(z) up to first order in z, or equivalently up to the term defining 
the /^-function for each model - remember that dA/dz = f3\ din H/dz. We may therefore 
distinguish models by plugging in the corresponding /3-function in this formula. 

In Fig. ||a the confidence region in the (0° , (1/ p c )dA/ ' dz | 2 =o) space for a flat universe 
is shown. It does by no means discard the running of the cosmological constant. In some 
sense this is similar to what happens with quintessence models where data are compatible 
with a slow evolution of a scalar field, but here we are testing the evolution of a density 
component without any variation in the equation of state. When we restrict ourselves to 
the constant cosmological term, dA/dz = 0, we recover results similar to the standard ones, 
namely tt° M ~ 0.3 ± 0.1, £l° A ~ 0.7 ± 0.1 for a flat geometry (Cf. Fig.ga). In this case the 
minimum % 2 value lies at £l° M = 0.29, Q A = 0.71. But in general we see from Fig. |8]a that 
this commonly accepted value may undergo a wide variation when the possible running of 
the cosmological constant is taken into account. 
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Figure 8: (a) Confidence regions for the current cosmological constant density and its variation 



with redshift for a flat universe- Cf. Eq. (6.10). From darkest to lightest the regions correspond to 
1(7, 2cr and 3cr respectively. A running CC clearly provides a wider interval of possible values for 
f2° with the same probability; (b) Confidence regions for an d the representative parameter v of 
our model with running CC. 



After this generic observation, we turn to our model and apply a x 2 test to obtain 
confidence regions in the (Jl^, ^)-plane for a flat universe (Fig.||b). We see there that 
)° 

A 



does not vary significantly with different values of v. This is because the sample of 



- 35 - 



supernovae used has redshift up to 0.83 (Cf. Fig.0) and the average is z ~ 0.5. However, 
as it was already observed in Fig. [l], the effect becomes more important as we go to higher 
z, and we will thus be able to draw more trustworthy conclusions when we analyze SNAP 
data, which will reach up to z = 1.7. 



7. SNAP and the running cosmological constant 

After the discovery of the accelerated expansion of the universe 0, [| a new satellite ob- 
servatory (SuperNova Acceleration Probe) was proposed to determine the nature of the 
dark energy cause of the acceleration. The SNAP collaboration aims to obtain spectra and 



photometry for 2,000 supernovae already in the first year of mission [34]. The distribution 
of supernovae will have a maximum in the interval 0.2 < z < 1.2 where according to the 
present observed rates 1800 supernovae should be found. A smaller number of data is 
expected to be obtained up to a redshift of 1.7 (see more details in |34|| ). 

We suppose that supernovae are equidistant within every redshift interval, and obtain 
the magnitude for these supernovae within the cosmological model presented in this work. 
In each case observational gaussian errors are added to these values taking into account 
the systematic uncertainties, which are estimated to be zero at z=0 and 0.02 at z=1.5 
following, for instance J55|]: <j sys = z (0.02/1.5) and the intrinsic dispersion of supernovae 
after the corresponding calibrations &i n t r = 0.15. 

After fitting the models to their simulated data we have an idea of the precision with 
which we will determine the parameters and the possibility of distinguishing among models. 
Moreover, when the first results of SNAP will become available, the value of Vt® M should 
already be known to a high precision. We can therefore adapt the present study to this 
foreseeable situation and use a Bayesian analysis in order to determine the probability of 
our parameters given some prior information. Details of the statistical methods used are 
shown in Appendix 2. 

7.1 High — z samples and the determination of v 

With the SNAP sample we will have, at high redshift, groups of 20 supernovae per bin 
of width 0.01 in redshift. This means that the uncertainty in the magnitude at each bin 
will be of the order of only 0.03 magnitudes for the present intrinsic dispersion of 0.15 
magnitudes. 

We can take one of this bins at relatively low redshift {z\ = 0.5) and another one at 
higher redshift (z% = 1.0 or 1.5) to look at the intersection between the allowed regions in 
the parameter space (Fig. 0). This intersection represents how the intrinsic dispersion in 
magnitude translates into the parameters according to (AJ7a x AS1 m) oc (cr^ x a!j%) (see 
for instance pq]). 

Figure |9] shows the results when the running of the cosmological constant is considered. 
In the left panel the bands of constant magnitude are plotted in the usual space of cosmo- 
logical parameters. The cosmological model used is one with = 0.3, n° A = 0.7, v = 2i/ 
and the bands give a similar intersection as in a standard model with VP M = 0.3, fi9 = 0.7 
(Fig. 2 in |56| ), so the present-day parameters will be determined under the same conditions 
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Figure 9: Bands of constant magnitude in the parameter space, where magnitudes are obtained for 
a universe with fl° M = 0.3, = 0.7, v — 2vq. The plot on the left shows the intersection between 
the bands coming from two supernovae at z = 0.5 and z — 1.0 on the (£1%, £1" )-plane. On the right 
the plot is on the z/)-plane where the two supernovae need to be more distant in redshift in 
order to equally restrict the confidence region. 



in both kinds of models. But we can also see how the cosmological index is determined 
by the SNAP measurements. In order to obtain a reduced confidence region in the (0° , u) 
plane, we need data more separated in redshift than in the previous case. This is be- 
cause the interesting models with standard (non- variable) cosmological constant have their 
maximum difference at z ~ 0.8 (this value is that obtained for the difference between a 
= 0.3, = 0.7 universe and a Q° M = 0.2, Q° A = universe). Models with a running 
cosmological constant differ from the ones without such running CC at higher redshift 



{z > 1) and the difference grows with z. Thus, as we see in Fig. 10, we need a sizeable 
group of supernovae at the highest redshift reachable by SNAP while a strong concentra- 
tion of data around z = 1 is needed to distinguish models with and without cosmological 
constant. 



7.2 Numerical simulations for SNAP and other high z SNe la samples 

We now use the distribution of data predicted by the SNAP collaboration and some al- 
ternative distributions to obtain the confidence regions for the representative parameter of 
the model, v. However, as we have done with the current data, we first use the distinctive 
characteristic of the variable CC models, namely the presence of a term dA/dz different 
from zero. The order of magnitude of these derivatives in our model model can be directly 
read from Fig. ||d. At the present epoch (1/ p®)dk/dz = 3v£l° M ~ 0.05 for v = 2vq. As 
seen in the figures, Eq. ( |6.10| ) is not a sufficiently good approximation in this case since the 
Taylor expansion is only significant at first order. Therefore, in our model we would have 
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Figure 10: Deviations of the distance modulus m — M with respect to the case of a standard 
cosmological constant with fij^ — 0-3, = 0-7 or equivalently v = 0. Two fully academic situations 
are shown: (fij^, f2 A ) = (0, 0) -Milne's open universe- and a pure de Sitter's universe (0, 0.5), both 
with no matter, in order to show the scope of the effects. Important differences appear for z > 0.5. 
The other curves in the graphic show the case of the more realistic model presented in the paper, 
for various values of v. Differences become important at z > 1. 

to use the second derivative as well in order to study the change of slope at high redshift. 
The other parameters used in the simulation are a flat universe with 09, = 0.3, A = 0.7. 

For such a universe we obtain, after the SNAP project, the results in Fig. [ll]. The left 
panel shows the confidence regions in the (0 A , (1/ ' p G c )dk/dz | 2 =o) plane when a flat universe 
is considered, and the right panel adds a gaussian prior on with aa M = 0.03. These 
restrictions in the analysis are justified by recent results of the WMAP experiment Q. 
Thus, with all these considerations, (l/p®)dA/dz at the present epoch would be obtained 
to within ±0.2 which could distinguish models like the one presented here from the standard 
cosmological constant model. 

The results for the parameter v are shown in Fig. Fits are made under the same 
assumptions as in the previous case and so the results represent the expected accuracy in the 
parameters using all the information to come. Thus we determine v to ±0.06 for v = 0.1. 
We must take into account, however, that the SNAP distribution of data is optimized in 
order to distinguish models that have their maximum difference at a redshift near 1. We 
can, therefore, try other distributions designed from the idea that the difference between 
models grows with redshift, as it indeed happens in our case. The first distribution is very 
similar to the SNAP one but with most of the data homogeneously distributed between 
z = 0.2 and z = 1.7. This second distribution extends data up to redshift z = 2. This is 
in part justified by data which are being obtained within the GOODS and HST Treasury 
Program |35| . 
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Figure 11: Confidence regions in the (£1%, (l/p®)dA/dz \ z =o) plane for a flat universe. The left 
(right) plot has been obtained without (with) restriction (i.e. prior) in the knowledge of Sl° M . Data 
used in the fits are those obtained with the SNAP distribution of 2,000 supernovae as explained in 
Section 0. 

The results from the various distributions are shown in Fig. ^ and quantified in Table)!]. 
The latter also compares the determination of the parameters representing an RG evolution 
of the cosmological constant (i.e. v ^ 0), and the evolution of the equation of state 
(see Eq. ( |2.9D ) for an alternative dark energy model. We see that the evolution of the 
cosmological constant can be determined with a relative error of 60% for a value of v = 0.1. 
Distributions of data other than SNAP (Distr.l and Distr.2 in Table p]) allow to reduce 
this error to 20% with v = 0.1, but determinations remain poor for small values of v, 
which on the other hand would be the natural ones. We are then in a situation similar to 
the determination of the evolution of the equation of state, Eq. Q2.9| ). Although wq can be 
obtained with a precision of 2% for wq ~ — 1, its evolution w\ is only determined with great 



uncertainty. For instance, the evolution of w for SUGRA described by [57| is expected to 
be determined with a relative error in w a ~ 2w\ of 43 %. In our case, as for other dark 
energy candidates, very higb-z samples as the ones being collected by GOODS and HST 
Treasury program will significantly contribute to the precise determination of v. 



8. Conclusions 

To summarize, we have constructed a semiclassical FLRW model with variable CC at the 
present cosmic scale. The variation of the vacuum energy is provided without introducing 
special scalar (quintessence-like) fields and is completely caused by quantum effects of 
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Figure 12: la and 2a confidence regions for the basic parameter v of our model and the current 
value of f^A- Left panels shows the results with SNAP data in a flat universe whereas right panels 
add a restriction (prior) on = 0.30 ± 0.03. The distinction is only made at la level within the 
range of v with physical meaning. 
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Figure 13: Confidence regions in the (£1^, v) space with data distributions (see Table [y) alternative 
to the standard SNAP one. This way provides a better determination of the parameters. 

vacuum. The evolution of the CC is due to the Renormalization Group running triggered 
by the smooth decoupling of the massive fields at low energies, while the RG scale fi being 
associated to the Hubble parameter, H, at the corresponding epoch. Although the f3\- 
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Distribution 


Data 


Prior au M 


e 


0-0 


SNAP 

(1 year) 


50 SNe < z < 0.2 
1800 SNe 0.2 < z < 1.2 
50 SNe 1.2 < z < 1.4 
15 SNe 1.4 < z < 1.7 


None 


v = 0.1 


±0.10 


SNAP 


as above 


0.03 


v = 0.1 


±0.06 


SNAP 


3 years 


None 


v = 0.1 


±0.06 


SNAP 


3 years 


0.03 


v = 0.1 


±0.04 


Distr.l 


50 blNe < z < 0.2 
2000 SNe 0.2 < z < 1.7 


0.03 


v = 0.1 


±0.05 


Distr.2 


250 SNe < z < 1 
1750 SNe 1< z < 2 


0.03 


v = 0.1 


±0.02 


SNAP 


3 years 


0.03 


w = -1 


±0.02 


SNAP 






w a = 0.58 


±0.25 



Table 1: Determination of the parameters with SNAP data and with other two distributions. In 
all cases the cosmology used for the analysis is a flat universe, and when a prior on is added we 
use a central value of CIm = 0.3. For comparison we show the parameters representing variations 
in the equation of state (for SUGRA in Q). 



function itself is proportional to the fourth power of the masses, the decoupling does still 
introduce an inverse power suppression by the heavy masses, and thus one is left with a 
residual quadratic law [3\ ~ H 2 M 2 . The effective scale M summarizes the presence of 
the heaviest degrees of freedom available. This peculiar form of decoupling stems directly 
from: i) the Appelquist-Carazzone (AC)-decoupling theorem, ii) general covariance, and 
also from iii) the non-fine-tuning hypothesis between the n = 1 ("soft decoupling") terms 
of /3a (Cf. Eq. (3.1)) whereby the overall coefficient of the quadratic contribution H 2 M 2 
does not vanish. This particular form of decoupling is a specific characteristic of the CC 
because it is of dimension four. There is no other parameter either in the SM or in the 
GUT models with such property. 

In constructing the cosmological model we have explored the possibility that the hea- 
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viest d.o.f. may be associated to particles having the masses just below the Planck scale. 
This assumption is essential to implement the soft decoupling hypothesis within the [i = H 
setting. Indeed, the present value of HqM 2 is just of the order of the CC, and therefore 
it insures a smooth running of the cosmological term around the present time. In this set- 
ting the /^-function has only one arbitrary parameter v ( 4.1Q ), and as a result the model 
has an essential predictive power. In general we expect \v\ <C 1 from phenomenological 
considerations, mainly based on the most conservative hypotheses on the nucleosynthesis 
framework. Furthermore, for \v\ <C 1, we insure the absence of the trans-Planckian ener- 
gies. It is not completely clear how much the cosmological index v can approach the value 
1 from below. Despite the values v < 1 are not completely ruled out, it is not clear that 
the conditions for the nucleosynthesis could be safe. However, if accepted, it opens up the 
possibility that the CC could be in transit from a A > regime into a future A < one, 
which would be a welcome feature for string/M-theory. 

Our RG model offers the possibility to explore the existence of the sub-Planck physics 
in direct cosmological experiments, such as SNAP (and the very high-z SNe la data to be 
obtained with HST). For example, for the flat FLRW case and a moderate and positive 
value v ~ 10 -2 , we predict an increase of 10 — 20% in the value of Q\ at redshifts z > 1 — 1.5 
perfectly reachable by SNAP. For similar, but negative, values of v we predict that the CC 
should become negative beyond redshifts z > 2. For v ~ 0.1 corrections to some FLRW 
cosmological parameters become as large as 50% or more. In general, this model has a wide 
spectrum of implications that could be tested by SNAP, even for fairly moderate values 
of v compatible with the most conservative bounds from nucleosynthesis. The simulations 
of the SNAP data in Section 7 show that the single parameter v of our model could be 
pinned down, for v < 0.1, with a precision of around 20 — 60%. Although this accuracy 
is not very high, the sign of the parameter could be determined and the effects would be 
manifest. Undoubtedly, this would be a good starting point to identify the presence of 
quantum corrections to the FLRW classical cosmology. 

The semiclassical FLRW model that we have proposed here explains the variation of 
the CC due to the "relic" quantum effects associated to the decoupling of the heaviest 
degrees of freedom below the Planck scale, and suggest that a time dependence of the CC 
may be achieved without resorting to scalar fields mimicking the cosmological term or to 
modifications of the structure of the SM of the strong and electroweak interactions and/or 
of the gravitational interactions. At the same time, our model provides a phenomenological 
parametrization for possible correlated deviations of the cosmological equations for matter 
density and dark energy. 

From a more fundamental point of view, we have shown that the cosmological con- 
stant and the matter density may evolve in a correlated way due to quantum effects without 
resorting to exceedingly exotic frameworks. Using the RG as a basic QFT tool and ex- 
trapolating the standard AC law of decoupling to the CC case, we have straightforwardly 
predicted a cubic redshift evolution law for the CC and the matter density. The larger the 
redshift that we can eventually explore the larger the effects that we predict. Most impor- 
tant, we have shown that this cubic law can be thoroughly tested by the next generation 
of cosmological measurements, which will be able to reach depths up to z ~ 2. 
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9. Appendix 1. On the non-local terms in the quantum corrections. 

In this appendix we shall go beyond the phenomenological presentation used throughout 
the main text and briefly discuss some theoretical aspects of our QFT framework. Despite 
the cosmological model developed in the main text of the paper looks appealing due to its 
relative simplicity and relation to the QFT corrections, it is not immediately clear which 
level of credibility can be granted to the hypothesis of the "soft" low-energy decoupling 
which we used here. As it was already remarked in the main text, the physical interpretation 
of the renormalization group in curved space-time requires the formalism beyond the limits 
of the well-known standard techniques. These techniques are essentially based on the 



minimal subtraction MS scheme of renormalization |43|, 44, [5^, 25 1, which does not admit 
to observe the decoupling. 

The physical interpretation of the renormalization group in the higher derivative sector 
can be achieved through the calculation of the polarization operator of gravitons arising 
from the matter loops in linearized gravity jh]]. In this case one can perform the cal- 
culations in the physical mass-dependent renormalization scheme (specifically in the mo- 
mentum subtraction scheme, in which \i is traded for an Euclidean momentum p), obtain 
explicit expressions for the physical beta-functions and observe the decoupling of massive 
particles at low energies. In the mass-dependent scheme one has direct control of the 
functional dependence of the /3-functions on the masses of all the fields, i.e. one knows 
explicitly the functions Pk( m i'-,P 2 /Mj) for the higher derivative terms in (|2.2[) . Quite 
reassuring is the fact that these momentum subtraction /3-functions boil down, in the UV 



limit, to the corresponding /3-functions in the MS scheme [46]. 

Unfortunately, as it was realized in |^ |, this program fails when we attempt to derive 
the physical (3- function for the CC (and also for the inverse Newton constant 1/G). In 
fact, within the momentum subtraction scheme we fail to see the /3-function for CC and 
for the inverse Newton constant 1/G at all, so there is no correspondence with the MS 
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scheme in these sectors. The origin of the problem is that the standard renormalization 
group running always corresponds to the non-local insertions in the effective action. For 
example, in case of QED the one-loop terms look, in the UV, like 

r QED = J V=~9 { - \f, v f^ + ^i^, log (-° ) F*> } . (9.i) 

In the MS-scheme of renormalization the total /x-independence of the effective action (see, 
e.g. p5fl ) is provided by the running of the coupling e = e(/x). This running is a physical 
effect which has been observed experimentally [59]. As noted before, if we wish to describe 
low energy effects p 2 <C m 2 , we need a more physical subtraction scheme where the role 
of fi is played by the momentum p. Then we expect decoupling terms of the form p 2 /m 2 , 
which manifest as additional □ / m 2 terms (without logarithms) in the effective action ([O]) . 
Similarly, in the higher derivative gravitational sector we have the terms of the form 

Thd = J d 4 xy/^ {C^/3 h(D, m 2 )C^ + Rk 2 (a,m 2 )R}, (9.2) 

where C^ va p is the Weyl tensor, R is a scalar curvature, m is any particle mass and k\ , ki 



are definite non-local functions [46] which in the high energy region include ln(D/m 2 ). As 
in the QED case, at low energy these functions contain the decoupling terms D/m 2 without 
logs. When the derivative operators act on the conformal factor a(t), they produce terms 
proportional to the powers of the Hubble parameter H and its derivatives. Since the 
derivative H of the Hubble parameter has the same order of magnitude as H 2 , for the 
sake of simplicity we can consider only the powers of H. It is clear that the finite order 
non-localities can be described, at low energies, by the local expansions with the tensor 
structure similar to the one of (|9.2| ). For instance, making the expansion in the powers of 
□/m 2 , we meet the quadratic law of the low-energy decoupling. 

Let us now consider the cosmological and Einstein-Hilbert terms. In the MS-scheme 
the /x-dependence of the form of the effective action is compensated by the running of the 
parameter A, similar to the QED case and the running e(p). As a result the overall effective 
action with a running CC is scale-independent, as it has to be. We are able, in principle, 
to observe the CC only at very small energies, and naturally expect to have much weaker 
running because of the low-energy decoupling. However, the quantitative investigation of 
the phenomena of decoupling for the CC met very serious difficulties, starting from the 



physical interpretation of the renormalization group in this case. As noticed in [46], the 
non-local insertions can not be done in the original gravitational terms, because acting on 
A by □ we get zero and □ acting on R gives a total derivative. The insertion can be 
done only into the non-local terms with similar global scaling, e.g. 

d 4 x^iv(^f) XT an d / d^x^Q R^(^)r^ (9.3) 

for the CC and inverse Newton constant respectively. For the effective action of the quan- 
tized massive fields these terms do not take place and hence the framework of linearized 
gravity is not appropriate for the physical interpretation of the renormalization group for 
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the CC and 1/G in curved space. Does it mean that there is no other framework where 
this interpretation should appear in a natural way? From our point of view it is not really 
so. Let us remind that: 

1. The renormalization group for the CC and 1/G can be formulated in a very robust 
way in curved space if we use the MS renormalization scheme. Therefore, we can expect 
the coincidence with the physical RG in an appropriate framework 9 . 

2. The linearized gravity approach is essentially based on the expansion near the flat 
background. But when quantizing massive fields we have to introduce the vacuum CC from 
the very beginning in order to provide renormalizability. In this case the flat Minkowski 
space-time is not a solution of the classical equations of motion and the whole approach 
based on the linearized gravity is formally illegal. Perhaps this is the source of the problem. 

In order to see the renormalization group and decoupling for the CC and 1/G one has 
to perform calculations of the polarization operator or vertices on the non-flat dynamical 
background. In this case there will be no relation between the covariant expansion in curva- 
tures and the gravitational perturbations around the given metric. From the point of view 
of the linearized gravity this means that we can not achieve the physical interpretation of 
the renormalization group via the gravitational perturbations around the flat metric back- 
ground, because all the effect is essentially non-perturbative. In some sense the situation 
is similar to the low-energy QCD, where one has to rely on the non-perturbative meth- 
ods or reformulate the perturbative calculations using, e.g. the renormalon technique. It 
might happen that the resummation of the perturbative series coming from both quantum 
loop expansion and the gravitational perturbations may lead to the visible renormalization 
group effect, but at present it is unclear how this could be checked. 

Let us consider again the linearized gravity approach. The non-localities of the effective 
action of the massive fields are related to the insertions of the "massive" Green functions 
1/(D + m 2 ) . Then, making the expansion in the powers of D/m 2 , we meet, typically, 
much stronger decoupling than the standard quadratic one. However, since we are dealing 
with an infinite multiple series, the correspondence with the MS scheme means that some 
resummation of the series in curvature produces the massless 1/D insertion which may 
eventually lead to the soft decoupling and to the UV correspondence with the MS scheme. 

At present there is no method of calculations on the non-flat background compatible 
with the physical renormalization scheme. In this situation our phenomenological approach 
looks rather justified. The covariance of the effective action forbids the first order in H 
corrections. Then, recalling that the UV contribution to j3\ from a particle of mass m 
is /?a ~ "i 4 _ Cf. Section H)- the low-energy contributions to /3a must be supressed by, at 
least, the factor 

U/m 2 ~ H 2 /m 2 , (9.4) 

and hence the overall low-energy /3a acquires the form H 2 m 2 , Eq. (|3^). The same form 
of decoupling can be expected for the induced counterpart Aj n d of the CC. The renor- 

9 Let us notice that the renormalization group for the CC and 1/G is nontrivial within the non- 
perturbative approach - see, e.g. [^o[. This supports our point of view on the scheme-dependence of 
the negative results of [Jm. 
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malization group equations for the vacuum and induced CC's are independent, and the 
relation emerges only at the moment when we choose the initial point of the RG trajectory 



for the vacuum counterpart [22, p3|] . In fact, the Feynman diagrams corresponding to the 
kind are > a t the cosmic scale, very similar to the vacuum ones, because the Higgs field 
(f> is in the vacuum state, < >^ 0, and the non-zero contribution of these diagrams is 
exclusively due to the momenta coming from the graviton external lines. Finally, there 
is no real need to distinguish induced and vacuum CC's in the present context. In both 
cases the absence of the n = order (non-supressed) contributions in (|3.1|) is required by 
the apparent correctness of the Einstein equations and the smallness of the observable CC. 
There is no guarantee that the IR supression does not take a much stronger form than the 
one (soft quadratic decoupling) assumed here. In the last case the future experiments and 
observations will demonstrate a perfectly constant CC 10 . The same should happen if the 
high-energy spectrum does not include the particles with the masses comparable to Mp. In 
both these cases the effect of running would be negligible. Then, since our renormalization- 
group based model gives very definite prediction for many observables (see Sections 4,5), it 
should be, sooner or later, distinguished from the standard model with the z-independent 
CC or maybe even from the quintessence models. 

One could wonder whether the non-local effects behind the renormalization group are 
compatible with the standard energy conservation equation which we used in the text. In 
fact, the covariant form of the conservation law < V a Tlf >= just reflects the covariance 
of the effective action and therefore does not depend on the non-localities which are always 
present in the quantum corrections. Furthermore, in a situation where the energy scale 
associated to the metric derivatives is very small, the effect of non-local terms can be, 
according to our model, presented in a compact form due to the renormalization group. 
In the leading order of magnitude one can use the renormalization-group induced terms 
to represent the leading non-local part of the effective action. In this case these terms 
should be many orders of magnitude larger than the ones corresponding to the next order 
approximation, and the conservation law should be valid especially for the renormalization 
group induced terms. 

Let us finally remark on the running in other sectors. For the inverse Newton constant 
1/G it is irrelevant because 1/G ~ Mp is very large and the effect of the running is 
relatively small as it was demonstrated in [j22jj. The relevance of the higher derivative 
terms at low energies is supposed to be negligible, as it was recently shown in [6C] for the 



perturbations of the conformal factor. Indeed, for the tensor perturbations the situation 
may be much more complicated (see, e.g. ffHJ and references therein). In contrast, the soft 
decoupling for the CC does matter because the CC is very small and any running, even 
the very small one that we get, is sufficient to produce some measurable effect at large 
redshifts. For precisely this reason we have to use the Friedmann equation, in combination 
with the full conservation law including the variable cosmological term, and no other terms. 



10 In case there is an additional source of the dark energy (e.g. quintessence), the variable dark energy 
will manifest itself. 



- 46 - 



10. Appendix 2. Bayesian analysis. 



In the Bayesian approach we can determine the probability density function for the cosmo- 
logical parameters (fijwj ^a) and model parameters {9 in general), conditional on the data 
(D). We know from Bayes's theorem that this quantity, p(£Im, ^a, 0\D), is obtained from 
the probability of the data conditional on the model and from the a priori information on 
the parameters and on the data: 

P (n M , n A , e\3) = p( ° |iiM ' ° A ' *V ( " M ' ■ (mi) 

p(D) 

Data coming from supenovae are effective magnitudes in the peak of the light curve, 
mf given in Eq. ( |6.8| ), and in this analysis 9 = (Ai, a,v). 

We assume that ^ are independent and each one follows a gaussian distribution. 
Then we obtain a probability density function for the data conditional on the model that 
is a product of gaussians. Hence the exponent is proportional to the sum that we have 
defined as x 2 i n Eq. Q6.9|), 



K^fl^M, fi A , 9) = I I exp [ p—L- I oc exp I ) . (10.2 



n M , U A , 9) = JJ * — exp [ 

i ,/2vra, 2 V 




As in the x 2 analysis we follow the results from [54] and we reduce 9 only to the 
parameter of the model (v in our case). Since we concentrate on the flat case we only 
introduce a gaussian prior on Qm- Using that p{m e /^) is constant too we obtain 

p(n „ ) = ^ exp (4 (1 „. 3) 



A I 



M 



P iO u SU.,\m?f) = lexp L£\ exp I , (i .4) 



where N is the normalization given by 



/>oo roo r-oo / 2 

N= dn M I dQ A dvexp f ) exp | _i ^l_Z^MI |. (i .5) 

JO J —oo J —oo V 
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